/ 



DOE/JPL-9 565 25-8 5/05 
9950-1181 


DRL No. 198 

DRD Line Item No. SE-8 


D0E/JPL-956525-85/8 
Distribution Category UC-63 


CAT. W 


//r- 


V37 ’3y 


HEAVY DOPING EFFECTS IN HIGH EFFICIENCY SILICON SOLAR CELLS 




Annual Report 


For Period Covering 
June 24, 1984 - June 23, 1985 


By: 

F. A. Lindholm and A. Neugroschel 
UNIVERSITY OF FLORIDA 
Department of Electrical Engineering 
Gainesville, FL 32611 


JPL Contract No. 956525 


The JPL Flat-Plate Solar Array Project is sponsored by the U. S. Department of 
energy and forms part of the Solar Photovoltaic Conversion Program to initiate 
a major effort toward the development of low-cost solar arrays. This work was 
performed for the Jet Propulsion Laboratory, California Institute of 
Technology, by agreement between NASA and DOE. 


(NASA-CR-179954) HEAVY DOTING EFFECTS IN 
HIGH EFFICIENCY SILICON SCLftE CEILS Annual 
fieport, 24 Jun 1964 - 23 Jun. 1965 (Florida 
Univ. 1 109 p CSCL 10, 


N6 7-1 3853 


Unclas 
G3/44 43631 


This work was performed for the Jet Propulsion 
Laboratory, California Institute of Technology, and was 
sponsored by the U. S. Department of Energy through an 
agreement with the National Aeronautic and Space 
Admi ni strati on . 

This report was prepared as an account of work 
sponsored by an agency of the United States Government. 
Neither the United States Government nor any agency thereof, 
nor any of their employees, makes any warranty, express or 
implied, or assumes any legal liability or responsibility 
for the accuracy, completeness, or usefulness of any 
information, apparatus, product, or process disclosed, or 
represents that its use would not infringe privately owned 
rights. 

Reference herein to any specific commercial product, 
process, or service by trade name, trademark, manufacturer, 
or otherwise, does not necessarily constitute or imply its 
endorsement, recommendation, or favoring by the United 
States Government or any agency thereof. The views and 
opinions of authors expressed herein do not necessarily 
state or reflect those of the United States Government or 
any agency thereof. 



TABLE OF CONTENTS 


SUMMARY (EXTENDED ABSTRACT) iv 

CHAPTER 

ONE INTRODUCTION 1-1 

TWO HEAVILY DOPED POLYSILICON-CONTACT SOLAR CELLS 2-1 

THREE EXPERIMENTAL STUDY OF THE MINORITY CARRIER TRANSPORT AT THE 

POLYSILICON-MONOSILICON INTERFACE 3-1 

FOUR SOLUTION OF THE CONTINUITY EQUATION IN PLANAR SYMMETRY CASES 

AND ASSESSMENT OF PHOTOLUMINESCENCE DECAY 4-1 

FIVE TRAP CONTROLLED MINORITY-CARRIER MOBILITY IN HEAVILY DOPED 

SILICON 5-1 

SIX DETAILED INTERPRETATION OF PHOTOCONDUCTIVITY DECAY 

RESPONSE FOR LIFETIME DETERMINATION 6-1 

SEVEN GENERALIZED RECIPROCITY THEOREM FOR SEMICONDUCTOR DEVICES 7-1 

EIGHT MINORITY-CARRIER DIFFUSION COEFFICIENTS AND MOBILITIES 

IN SILICON 8-1 

NINE DISCUSSION AND RECOMMENDATION 9-1 


PME BUW NOT FHJ*D 


jii 


SUMMARY (EXTENDED ABSTRACT) 

We report the first use of a (silicon)/(heavily doped polysilicon)/ 

(metal) structure to replace the conventional high-low junction (or back- 
surface-field, BSF) structure of silicon solar cells. Compared with BSF and 
back-ohmic-contact control samples, the polysilicon-back solar cells show 
improvements in red spectral response and open-circuit voltage. Measurement 
reveals that a decrease in effective surface recombination velocity S is 
responsible for this improvement. Decreased S results for n-type (Si:As) 
polysilicon, consistent with past findings for bipolar transistors, and for 
p-type (Si : B) polysilicon, reported here for the first time. Though the 
present polysilicon-back solar cells are far from optimal, the results suggest 
a new class of designs for high efficiency silicon solar cells. Detailed 
technical reasons are advanced to support this view. 

We present the results of an experimental study designed to explore both 
qualitatively and quantitatively the mechanism of the improved current gain in 
bipolar transistors with polysilicon emitter contact. Polysilicon contacts 
were deposited and annealed at different conditions. The electrical 
properties were measured using p/n junction test structures that are much more 
sensitive to the contact properties than are bipolar transistors. A simple 
phenomenological model was used to correlate the structural properties with 
electrical measurements. Possible transport mechanisms are examined and 
conjectures are made about' upper bounds on transport parameters in the 
principle regions of the devices. The main conclusion of this study is that 
the minority carrier transport in the polycrystalline silicon is dominated by 
a highly disordered layer at the polysilicon-monosilicon interface 
characterized by very low minority-carrier mobility. The effective 
recombination velocity at the n + polysilicon - n + monosilicon interface was 
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found to be a strong function of fabrication conditions. The results indicate 
that the recombination velocity can be much smaller than 10^ cm/sec. 

The decay of excess carriers induced in a semiconductor sample by a laser 
pulse is studied. The decay rate is proved to be independent of the spatial 
dependence of the hole-electron generation rate provoked by the laser beam. 

The effect of surface recombination can be minimized by waiting until the 
response drops to about one-percent of its initial value. The accuracy of 
photoluminescence decay measurements of Auger lifetime is examined in view of 
this and other criteria developed here. 

The activation behavior of the minority-carrier mobility and diffusivity 
20 -3 

in heavily doped (~ 10 cm ) Si (Si :As) is investigated in the temperature 
range, 20 K to 350 K. Experimental results indicate that hole transitions 
between the valence band and localized shallow states give rise to the 
observed behavior. The activation energy is about 10 meV, which suggests that 
the localized states originate from band tails but does not rule out trapping 
at boron atoms in the compensated n + region. 

Photoconducti ve decay provides a method for estimating recombination 
lifetime in semiconductors. It often exhibits a single characteristic time. 
For samples having large recombination-center densities, previous experiments 
indicate that this characteristic time can considerably exceed the steady- 
state lifetime. We present apparently the first explanation of these 
results. The explanation involves the effect on the overall photoconducti ve 
decay of the multiple characteri Stic times that describe trapping and 
recombination through bound states for variations of the electrochemical 
potentials that are small relative to the thermal voltage. 

A reciprocity theorem is presented that relates the short circuit current 
of a device, induced by a carrier generation source, to the minority carrier 
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Fermi level in the dark. The basic relation is general under low injection. 

It holds for three dimensional devices with position dependent parameters 
(energy gap, electron affinity, mobility, etc.), and for transient or steady- 
state conditions. This theorem allows calculation of the internal quantum 
efficiency of a solar cell by using the analysis of the device in the dark. 
Other applications could involve measurements of various device parameters, 
interfacial surface recombination velocity at a polycrystalline silicon 
emitter contact, for example, by using steady-state or transient photon or 
mass-particle radiation. 

A new method for accurate measurement of minority-carrier diffusion 
coefficients in silicon is described. The method is based on a direct 
measurement of the minority-carrier transit time through a narrow region of 
the p-n junction diode. The minority-carrier mobility is obtained from the 
diffusion coefficient using the Einstein relation. The method is demonstrated 
on low-doped n- and -p-type Si (dopings ~ 10^ cm"^) and is compared with the 
literature data for the majority-carrier mobilities. The results show that in 
low-doped Si the electron minority- and -majority-carrier mobilities are 
comparable, but the hole minority-carrier mobility is significantly higher 
(~ 30%) than the corresponding majority-carrier value. The results confirm 
earlier data of Dziewior and Silber. 
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CHAPTER ONE 
INTRODUCTION 

This annual report describes technical findings resulting from JPL 
contract 956525 for the period June 24, 1984 to June 23, 1985. 

The purpose of this contract is to assess the status of critical 
parameters describing heavily doped silicon, to integrate this assessment into 
an effort aiming toward determination of these parameters and their 
incorporation in silicon solar cell design, and to develop associated 
measurement methods. The interest in heavily doped silicon extends to heavily 
doped polysilicon, particularly to its use as a quasi -blocking contact at 
surfaces. It extends also to measurement methods designed to explore dilutely 
doped silicon, which typically adjoins heavily doped layers in silicon solar 
cells. 

Professors Fred A. Lindholm and Arnost Neugroschel of the Electrical 
Engineering Department of the University of Florida serve as principal 
investigators for this contract. Several graduate students have participated 
significantly in the conduct of the research: T. W. Jung, J. J. Liou, 

K. Misiakos, and J. S. Park. Other collaborators have included Peter lies, 

M. Arienzo, B. Y. Hwang, Y. Komen, R. Issacson, M. Spitzer and Professors P. 

T. Landsberg and C. T. Sah. 

In the brief description to follow presently concerning the content of 
the chapters of this report, some of these people named above are indicated, 
through the appearance of their initials, as principal contributors to the 
chapter being described. 

Prior to outlining the content of each chapter, we list immediately below 
the journal papers (two invited) and the conference papers (including one 
invited plenary lecture) that have resulted from this contract research. Only 
those papers for the present reporting period, June 1984 to June 1985, are 
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listed. We include both papers published and papers accepted during this 
period. 




PAPERS SUPPORTED BY CONTRACT #956525 

F. A. Lindholm and J. J. Liou, "Improved determination of lifetime and surface 
recombination velocity by observing transient voltages," 18th IEEE 
Photovoltaic Specialists Conference, October 1985 (invited plenary lecture to 
be presented). ~~ ~~ ~~~ 

F. A. Lindholm, A. Neugroschel , M. Arienzo, and P. A. lies, "Front and back 
polysilicon-contacted solar cells," 18th IEEE Photovoltaic Specialists 
Conference , October 1985. 

K. Misiakos and F. A. Lindholm, "Generalized reciprocity theorem for 
semiconductor devices," J. Appl . Phys., tentatively scheduled for October 1, 
1985. 

K. Misiakos and F. A. Lindholm, "Toward a systematic design theory for silicon 
solar cells using optimization techniques". Solar Cells, Oct. 1985 (invited 
paper). 

A. Neugroschel, F. A. Lindholm, and C. T. Sah, "Trap controlled minority- 
carrier mobility in heavily doped silicon," Solar Cells, to appear in 
Sept. 1985. 

A. Neugroschel, "Minority-carrier diffusion coefficients and mobilities in 
silicon," IEEE Electron Device Letters , vol . EDL-6, pp. 425-427, August 1985. 

K. Misiakos, F. A. Lindholm, and A. Neugroschel, "Solution of the continuity 
equation in planar symetry cases and assessment of photoluminescence decay," 

J. Applied Physics , vol. 58, pp. 1647-1650, August 1985. 

F. A. Lindholm, A. Neugroschel, M. Arienzo, and P. A. lies, "Heavily doped 
polysilicon contact solar cells," IEEE Electron Device Letters , vol. EDL-6, 
pp. 363-365, July 1985. 

P. T. Landsberg, A. Neugroschel, F. A. Lindholm, and C. T. Sah, "A model for 
band-gap shrinkage in semiconductors with application to silicon," Physica 
Status Solidi (a) , to be published July 1985. 

F. A. Lindholm, "Uncertainties about the physical electronics of n + and p + 
silicon, with applications for solar cells". Solar Cells , vol. 12, pp. 131— 
140, June-July 1984 (invited paper). 

B. Y. Hwang and F. A. Lindholm, "Detailed interpretation of photoconductivity 
decay response for lifetime determination," Solar Cells , vol. 14, No. 2, 

pp. 187-190, May 1985. 
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A. Neugroschel , M. Arienzo, Y. Komen, and R. D. Isaac, "Experimental study of 
the minority carrier transport at the polysilicon monosilicon interface," IEEE 
Trans. Electron Devices , vol . ED-32, pp. 807-816, April. 

0. von Roos and F. A. Lindholm, "Steady state currents in semiconductor 
filaments or grains in case of large surface recombination velocity on lateral 
surfaces," J. Appl . Phys. , vol. 57, pp. 415-417, January 1985. 

This report proceeds as follows. Chapter Two (FL, AN, MA, PI) describes 
the first use of polycrystalline silicon as a passivant for the surfaces of 
silicon solar cells. Chapter Three (AN, MA, YK, RI) reports experimental 
findings about the mechanisms underlying this passivation — that is, this 
reduction of surface recombination velocity. For doping concentrations of 
about 10^/cm^, we inferred that the value of this velocity could be as small 
as 15 cm/s for small -area devices. Chapter Four (KM, FL, AN) provides a 
detailed mathematical assessment of the Auger-impact rate coefficients 
determined by the photoluminescent decay method (Dziewior and Schmid, 1977). 

We conclude that in principle the method can yield accurate values despite the 
presence of large surface recombination velocity and despite areal 
inhomogeneity of the incident laser beam. Chapter Five (AN, FL, CS) examines 
experimentally the credibility of the model for values of the minority-carrier 
mobility and diffusivity that are much smaller than the values of the majority 
carrier counterparts, at least for doping concentrations substantively above 
10^/cm^. This model involved trapping of minority carrier at bandtail or 
other shallow bound states associated with the minority-carrier band 
(Neugroschel and Lindholm, 1983). The results support this model. 

Chapter Six (BH, FL) provides a detailed interpretation of the photoconducti ve 
decay method for examining volume or bulk lifetime in silicon wafers. 

Chapter Seven (KM, FL) puts forward a generalized reciprocity theorem for 
semiconductor devices. This theorem allows calculation of the internal 
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quantum efficiency of a solar cell by using the analysis of the device in the 
dark. Chapter Eight (AN) sets forth experimental findings concerning the 
measurements of minority-carrier mobility and diffusion coefficient in 
dilutely doped silicon (about 10^/cm^). The results show that the hole 
minority-carrier mobility is 30% higher than the majority-carrier value. 
Design considerations immediately follow from these findings. Chapter Nine 
briefly discusses the results of this annual report. 
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CHAPTER TWO 

HEAVILY DOPED POLYSILICON-CONTACT SOLAR CELLS 
INTRODUCTION 

Highly conducting polysilicon contacts to crystalline silicon are widely 
used in bipolar integrated circuits [1]. Among other properties, these 
heterojunctions reduce recombination current by forming quasi-blocking 
contacts [2]. 

Our purpose is to describe experimental results and their implications 
deriving from the first use of such contacts in silicon solar cells. 
Specifically, in this paper we report the results of introducing such contacts 
at the back surface. The resulting solar cells show large improvements in red 
spectral response RSR and moderate improvements in open-circuit voltage V QC 
when compared with control cells made with the conventional back-surface-field 
(BSF) structure or with back-ohmic contact (BOC) structures. To assess the 
origin of these improvements, we measure the effective recombination velocity 
S of minority carriers in the monocrystalline silicon as they enter the 
polysilicon layer. This determines that S( pol ysi 1 icon back) < S(BSF 
controls), and we interpret this lowering of the recombination velocity to be 
the mechanism responsible for the improvement in V oc and RSR. The solar cells 
and controls described here have n-type and p-type base regions; the decrease 
in S resulting from p-type polysilicon (Si : B) is reported here for the first 
time. 


FABRICATION 

The substrate resistivity was 2 - 0.5 n-cm for both the p-type and n-type 
samples. The p-type substrates were chemically polished on the back side; the 
n-type substrates were polished chemically-mechanically. About 2000 A thick 
layers of n + and p + polysilicon was deposited (in-situ doped) at 670°C in an 
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atomospheric pressure CVD reactor. The n + -emitter (Si : P) was diffused at 
87S°C for about 12 min., and the p + -emitter (Si rB) was diffused at 9Ul)°C for 
30 min. During the polysilicon deposition and the emitter diffusion, CV0-Si02 
protected the other side of the sample. Out-diffusion from the doped 
polysilicon into the base occuring during the emitter diffusion creates a 
crystalline low-high junction of thickness <200 A [2]. The thinness of the 
resulting n + (or p + ) layer at the back suggests that negligible recombination 
occurs in its volume. 

The conventional Al-paste alloying technique formed the BSF region to a 
depth of about 3 to 5 pm for the p-base BSF controls. Phosphorous diffusion 
formed a one-micron thick BSF region for the n-base BSF controls. These 
n-base controls and back-ohmic-contact controls had an ohmic contact (Ti + Pd 
+ Ag) deposited and then annealed on the back surface. The p-base samples had 
a thin layer of A1 under the Ti + Pd + Ag, and the heating was below the 
Al-alloy temperature (~ 4UU° C) . 

MEASUREMENTS AND EXPERIMENTAL DATA 

We made measurements to compare the merit of the polysilicon-back-contact 
structure to that of the BSF and BOC controls. Our main object was to explore 
the influence of the back surface recombination velocity, which enters the 
boundary-value problem describing the solar cells. 

To obtain the boundary-value problem, we follow Shockley in thinking of 
the device as comprising two quasi -neutral regions separated by a junction 
space-charge region [3]. We focus on recombination in the volume and at the 
surface of tne quasi-neutral base. If we define the surface as the interface 
between the silicon and the polysilicon contact, characterized by S, we get a 
simple model in which the polysilicon recombination and transport mechanisms 
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are imbedded in S. For a BSF cell, we define the surface as that plane where 
the low-high (p/p+ or n/n+) crystalline junction begins. Then S characterizes 
recombination occuring in the transition and heavily doped regions. (The 
usual carrier density and voltage relations together with the continuity 
equation for minority carriers, completes the boundary-value problem.) 

Using an AMO spectrum, we measured the solar-cell parameters. Of these 
the one most closely linked to S is the red spectral response RSR, 

{\ > 6000 A). The data appear in Table 1. For a p-type base , the cells 
having a polysilicon contact show significant increases in V QC and RSR over 
those having an ohmic back contact or an aluminum-alloy back surface field. 
Note, however, that the alloy BSF controls showed no improvement compared with 
the BOC controls in V QC and only a slight increase in RSR. This indicates 
that the alloy process failed to yield an effective back-surface field. It is 
common experience that the alloy process does not consistently yield good BSF 
regions for p as low as 2 Q-cm, although it is reliable for higher base, 
resistivities. For the cells having an n-type base, both the polysilicon 
devices and the BSF controls, formed by P diffusion, showed significant and 
similar increases in RSR when compared with the response of the BOC 
controls . 

To assess whether one can attribute these results to smaller S, we 
measured S and diffusion length L using an MOS-switch version of the method of 
Ref. [4J. In these experiments, the silicon-polysilicon junction provides a 
blocking contact that is markedly superior to the low-high junction of the BSF 
control s (Table 1) . 
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DISCUSSION AND EXTENSIONS 


The most important columns in Table 1 pertain to the red spectral 
response and to S, which emphasize the improvements deriving from the 
polysilicon/silicon back contact. The improvement in V QC is less striking, 
partly because V QC depends not only on S but also on the collective effects of 
L, shunt conductance and recombination current in the junction space-charge 
region and in the quasi -neutral emitter. 

The polysilicon/silicon junction is also superior in other respects. 
Process control for BSF devices is critical for good performance. This is 
particularly true for the A1 alloy process commonly used to fabricate a p-type 
BSF region. As noted, this process does not reliably form good BSF regions 
for p ~ 2 Q-cm. For BSF regions formed by solid-state diffusion, the high 
temperatures required may degrade lifetime and the low-high junction may be 
shifted by subsequent thermal processing such as that used to form the p/n 
junction. Both diffusion and alloy processing consumes part of the substrate, 
which may pose difficulties for thin solar cells. In contrast the 
polysilicon-back structure derives from a non-critical deposition process at 
low temperature [2]; the post-deposition heat treatment also is not 
critical [2]. Polysilicon back processing consumes essentially none of the 
substrate. It can be controlled precisely and is reliable, as employment in 
bipolar integrated-circuit technology has shown. 

Consider now extensions of the present design. First, consider low-p 
substrates (0.5 Q-cm < p < 0.1 Q-cm, corresponding to doping concentrations 
~ 1017 cm” 3) , which reduce base-region volume recombination, thus tending to 
increase V QC and efficiency n» However, effective BSF regions on low-p 
substrates have never been made. In BSF solar cells, Sp * p - ^ for low- 
injection conditions [5J and for p > 0.5 Q-cm approximately. For 
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p < 0.5 n-cm,the potential barrier V(BSF) between the low- (~10^ cnT^) and 
high-doped sides of the BSF junction becomes ineffective: V(BSF) =* 
(kT/e)log[N(low)/N e f.p] < 2.3(kT/e) where N(low) is the doping concentration on 
the low-doped side and N e ff ~ cm"^ is the effective doping concentration 
on the high-doped side, corrected for band-gap shrinkage AEq. Thus from 
Ref. 5, S = (D/L) [N(low)/N e ^(high)] D/L = S(max), where the minority- 
carrier diffusivity D and diffusion length L pertain to the highly doped 
side. For a polysilicon-contacted solar cell, recent experiments [2] show 
that S max (poly junction) « S max (BSF). The mechanisms underlying the 
relatively low value of S max (poly junction) are still debated in the 
literature. The polysilicon of Ref. 2 and that for our devices of Table 1 was 
deposited without an intentional interfacial oxide, and Ref. 2 advances the 
view that the low value of S(max) derives from extremely low values of 0 in 
the highly disordered interfacial region. 

Thus we anticipate that polysilicon deposition can give low-S contacts on 
base regions having low resistivities (~ 0.1 ft-cm). From the results of 
Ref. 2, we anticipate S < 1000 cm/s. 

Attaining high efficiency (n) depends on reducing S at the front surface 
as well. Thermal oxides have proved effective over the surface, exclusive of 
that portion where ohmic contacts exist. Even though ohmic contacts 
constitute only 5% or less of the front surface, one should note that the 
increase in emitter recombination current caused by areal inhomogeneity 
greatly exceeds that calculated by consideration of the areal ratio of low-S 
to high-S regions [6,7]. To avoid the resultant drop in n. Green et al . [8] 
have used thin tunneling oxides between the silicon and the metal. 

Substitution of heavily doped polysilicon for the tunneling oxides may provide 
a better contacting structure with higher yield and longer life. 
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Sah has recently proposed detailed practical designs for controlling the 
recombination losses at the front, or emitter, surface [9]-[10J. He has 
proposed using thin and doped polysilicon layers between the metal and the 
heavily doped silicon surface of the emitter while using a thermal oxide over 
the emitter area not shadowed by metal. Employing the data of Ref. 2, he has 
calculated a practical maximum efficiency of above 23% provided by the low 
recombination losses at the front and back surfaces and assumed volume 
recombination only by Auger and radiative processes. 

Thus, we suggest elements of the following new class of silicon solar 
cell designs: 1) Use thin (=* 150 |im) base regions of resistivity ~ 0.1 Q-cm 

contacted on the back by heavily doped polysilicon; 2) Interleave the front 
surface with oxide and a polysilicon/metal system; 3) Use either an n-type or 
p-type silicon for the base region. 

As suggested by the work of Sah quoted above, we expect that these solar 
cells could have very high conversion efficiencies. The present work, 
together with the somewhat related efforts using SIPOS junctions to reduce 
surface recombination [11], constitute recent significant departures from 
conventional silicon solar cell designs. 
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TABLE 1 


Summary of measured parameters at 28° C, (no AR coating). The cells have 
area A = 4 cm^ and thickness of about 210 ± 10 pm except for cell 2N which was 
330 pm thick. The results are the average values from 10-20 cells. 


Cell 

Base 

type 

Back 

contact 

V 

(mV) 

I$Cr( X> 0.6pm) 
(mA) 

L 

(pm) 

S 

(cm/s) 

IP 

p 

ohmic 

573 

63 

310-350 

A A 

2P 

p 

BSF 

574 

64 

310-350 

4.2x10-5x10** 

5P 

P 

poly-Si 

BSF 

583 

67.7 

310-350 

1100-1500 

IN 

n 

ohmic 

566 

62.7 

190-250 


2N 

n 

BSF 

565 

65.2 

190-250 

700-1000 

5N 

n 

poly-Si 

591 

65.5 

190-250 

100-160 


BSF 
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CHAPTER THREE 

EXPERIMENTAL STUDY OF THE MINORITY CARRIER TRANSPORT AT 
THE POLYSILICON-MONOSILICON INTERFACE 

I. INTRODUCTION 

It was demonstrated [1-3] that the common-emitter current gain of bipolar 
transistors can be improved by replacing the emitter metal contact by a highly 
doped polysilicon layer. A number of works, both theoretical and 
experimental, investigating this effect, were published [1-8]. In all these 
studies, the analysis was done on the actual emitter structure of the bipolar 
transistor consisting of a thin (~ 2000 A) heavily-doped monocrystalline 
emitter region contacted by a heavily-doped polysilicon layer about 2000 A 
thick. The analysis of such a structure is difficult due to uncertainties in 
the parameters of both the monosilicon and polysilicon bulks and a poor 
understanding of the interface between the polysilicon and the single crystal 
silicon layer. These uncertainties are further augmented by the presence of 
the grain boundaries in the polysilicon. Thus, numerous assumptions are 
required to establish a minority-carrier transport model in the above system 
[4-8]. 

The region which is most difficult to characterize is the polysilicon- 
monosilicon interface. This region is very sensitive to the fabrication 
conditions and to the surface treatment before the CVD process of the 
polysilicon deposition [2,7]. A peak in the arsenic and phosphorus 
concentration was observed [7,9] as a result of the donor dopant segregation 
at the interface. The interface may also act as a natural sink for deep-level 
impurities and cause precipitation of interstitial oxygen from the bulk. The 
properties of the interface are very strongly influenced by a presence of a 
very thin insulating interfacial layer created by oxidizing silicon before the 
polysilicion deposition. This can be achieved by thermal oxidation [2] or by 
a chemical treatment [7]. 
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These and other variables make the evaluation of the relative importance 
of the polysilicon bulk and the interface very difficult, both theoretically 
and experimentally. De Graaf and de Groot C2] investigated the current gain 
of the n + /p/n transistors with an intentionally grown oxide layer at the 
polysilicon-monosilicon interface. They suggested that the current gain 
enhancement is controlled by the tunneling mechanism through the oxide 
layer. Ning and Isaac [3] have shown that, for the n + /p/n transistors with a 
polysilicon emitter contact with no intentional oxide layer at the 
polysilicon-monosilicon interface, the minority hole base current becomes a 
function of polysilicon thickness if the film thickness is less than about 
500 A. This indicates that the current gain enhancement is determined by the 
hole transport in the polysilcon film within about 500 A of the interface. 

The devices used in [3] with different thicknesses of the polysilicon film had 
identical polysilicon-monosilicon interfaces. 

In the present work we have examined the roles of the polysil icon- 
monosilicon interface and of the polysilicon bulk in much more detail than has 
been done previously. We intentionally emphasize the interface using devices 
with identical interface properties, while varying drastically the polysilicon 
bulk properties. This is achieved by using a bilayer structure of polysilicon 
films £9], described in detail in the next section. Our approach is, thus, 
different from that taken by other researchers [2,7J who varied the interface 
treatment while depositing the same polysilicon layer. 

The microstructure of the polysilicon film and of the interface was 
examined by transmission electron microscopy (TEM) and the distribution of 
arsenic in the polysilicon film and at the polysilicon-monosilicon interface 
was studied by means of secondary ion mass spectrometry (SIMS) [9], These 
structural properties of the polysilicon film and the interface are 
correlated, for the first time, with the electrical measurements. 
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To study the polysilicon-monosilicon interface we use a p/n junction test 
structure, rather than a bipolar transistor. The test structure and its 
fabrication are described in Section II. As is discussed later, this 
structure has a number of important advantages in comparison with the bipolar 
transistor and will be instrumental in developing the phenomenological 
transport model and revealing the dominant recombination mechanisms in the 
polysilicon. Sections III and IV. 

II. DESCRIPTION OF DEVICES AND FABRICATION 

The test structure, shown in Fig. 1, consists of a 0.06 Qcm p-type 

substrate with a thin (~ 0.8 pm) n-type epitaxial layer with doping density 
16 3 

Ngp - 10 cm grown on top. One kind of devices had an n + contact ring 
implanted; the second kind of devices on the same wafer does not have the 
ring. The purpose of the n + ring, which covers less than about 6% of the 
device area, is to assure a good ohmic contact with the relatively low-doped 
epitaxial layer. One half of each wafer was contacted by a 1500 A or 2500 A 
thick polysilicon film. Fig. 1(a). The deposition and the heat treatment 
after the deposition were done very carefully and were specifically designed 
for the purpose of enabling the separation of the role of the interface and 
the bulk on the current transport. Wafers were etched in a buffered HF 
solution prior to polysilicon deposition. It is very important to emphasize 
here that each kind of polysilicon device had a reference device with a metal 
A1 contact (both with and without the n + ring), Fig. 1(b), made on the other 
half of each wafer. The A1 contact was also deposited on top of the 
polysilicon film. The devices are separated from each other by an oxide layer 
around the periphery that minimizes the surface leakage, thus assuring that 
the measured current originates mostly from the bulk. The device areas ranged 
from 6 x 10"^ cm 2 to 1.2 x 10" 2 cm 2 . 
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The polysilicon films, deposited using an atmospheric pressure CVD 
reactor at 670°C, were one of the following: 

(a) 1500 A undoped polysilicon deposited without any external arsenic 
source (grain size about 100 A). 

(b) 1600 A in-situ arsenic doped (grain size larger than about 400 A). 

(c) Bilayer: 1500 A of undoped film followed by deposition of 1000 A 
in-situ doped film (grain size larger than about 300 A). 

Each of these three groups of devices was subjected to various post-deposition 

treatments. The deposition and annealing conditions are summarized in 

Table I. _ _ 

The arsenic-concentration profiles measured by SIMS [9] are shown in 

Fig. 2 for some of the devices from Table I. These profiles depict the 

average concentration of the cross section consisting of the grain-boundary 

region and the bulk. There is evidence [10] that in the in-situ doped films, 

about 30-50% of the arsenic dopant resides in the grain boundaries. Since the 

grain boundaries have smaller cross-sectional area than the adjoining grains, 

the arsenic concentration in the grain boundaries is larger than in the 

bulk. Based on the diffusion coefficients of arsenic in the polysilicon [9], 

we infer from the SIMS profiles of the in-situ doped films (Fig. 2.(a)), that 

20 -3 . . 21 

the arsenic concentration is about 5 x 10 cm in the grams and about 10 
cnT^ in the grain boundaries. The dopant in the as-deposited polysilicon film 
is largely inactive and is subsequently activated during the heat treatment 
step. The activation ratio (carrier concentration/dopant concentration) 
depends strongly on the annealing temperature and time [11,12]. For the 
in-situ doped sample 2BE annealed at 1000°C, 15 min, for example, the 
activation ratio is about 0.1 [11]. The activation ratio for the devices 
annealed at 900°C (devices 1C, 10) and at 800°C (device IB) is not exactly 
known, but it can be expected to be about 0.01 [11,12]. 
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For the wafers with the bi layer polysilicon film we have employed arsenic 

diffusion at low temperatures to obtain a preferential grain boundary 

diffusion from the top in-situ doped layer towards the interface, while 

minimizing the bulk diffusion [9]. The diffusivity of arsenic in the grain 

boundaries at T = 800 - 850°C is about three orders of magnitude greater than 

the diffusivity in the bulk. The SIMS profiles for some of the bilayer 

samples, before and after the heat treatments, are shown in Fig. 2(b). The 

arsenic concentration in the originally undoped layer depends upon the heat 

treatment. As an example, in the bilayer film of the sample 1J after heat 

treatment at 800°C for 64 hours, arsenic diffused through the grain boundaries 

and reached the interface [9]. From the SIMS profile, the average arsenic 

concentration in the grain boundaries is estimated to be about 5 x 

10 21 * cm"'*. However, in the sample 2BA with a heat treatment at 751)° C for 

8 hours the arsenic does not reach the polysilicon-monosilicon interface and 

the arsenic concentration drops from about 10 2 * cm"'* in tne in-situ doped part 

17 3 1 H 3 

of the film to the background concentration of about 10 cm - 10 cm 
[9]. This background level is due to the residual impurities in the 
deposition system and was detected in the undoped film. 

Although during the heat treatment arsenic diffuses mainly through the 
grain boundaries, some diffusion occurs from the grain boundaries to the bulk 
of the grains. To minimize this outdiffusion, the annealing conditions were 
tailored in such a way as to allow arsenic to just reach the interface. The 
outdiffused region is very narrow especially at the interface, where arsenic 
outdiffuses only during a fraction of the total heat treatment time. For 
example, for 800°C, 10 hrs, the arsenic concentration should drop from 2 x 
10 20 cm"'* to 10 17 cm""* in about only 100 A. During a fraction of the heat 
treatment time, arsenic will also diffuse into the epitaxial layer. The width 
of this diffusion should also be in the order of ~ 100 A. As a result of the 
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preferential grain boundary diffusion, a heavily doped thin interfacial region 
AVI is created at the metallurgical polysilicon-monosilicon interface. 

Figure 3 illustrates schematically the expected arsenic concentration profiles 
at the interface corresponding to the SIMS data data for the sample 1J from 
Fig. 2(b). 

Comparison of the SIMS profiles, Figs. 2(a) and 2(b), clearly shows that 
the in-situ doped and the bilayer devices have similar arsenic, and thus 
electron, concentrations within about 100-300 A thick interfacial layer AVI, 
but their bulk concentrations are different (see Fig. 3). Due to the 
activation ratio, the in-situ doped devices have a bulk electron concentration 

lO IQ o lO O 

of about N = 5 x 10 - 5 x 10 cm --5 in comparison with N « 10 cm for 

the bi layer samples in the originally undoped layer adjacent to the interface. 
As will be discussed in detail later, separation of the interfacial layer from 
the rest of the polysilicon film is also necessitated by very different 
microstructures in these two regions. 

Based on this similarity of the polysilicon-monosilicon interfacial 
layer, if the interfacial layer dominates the minority carrier transport, we 
expect the in-situ and the bilayer samples exhibit similar electrical 
properties. If, however, the bulk properties of the polysilicon determine the 
carrier transport, then we would expect the in-situ doped and the bilayer 
devices to have different electrical properties. 

III. MODEL DEVELOPMENT 

In this section we develop a simple phenomenological model for the 
minority-carrier transport with the aim of identifying the dominant 
recombination mechanism in the polysilicon-monosilicon structure. 
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Figure 4(a) shows the test structure under forward bias. For the present 
purpose we consider the polysilicon-monosilicon interface (Fig. 3) to be a 
thin layer of width AW ~ 100 A that includes the high-low junction in the 
epitaxial layer, grain-boundary space-charge-region (SCR) at the metallurgical 
interface and the highly doped region adjacent to the metallurgical interface 
in the polysilicon. We assume that this region can be characterized by a 
single effective time constant -t p (int). This assumption simplifies the 
analysis and will yield useful results. 

Using a charge-control representation for the hole current J p and the 
electron current J n at x=0, we write 


J p (0) 


Qp(epi) Q p (int) Q p (poly) 

-CpIepT ) + T p (int) + T p (poly) + J p (W E ) 


J n (0) 


Q n (sub) 

'tp(sub) 


( 1 ) 

( 2 ) 


J = J p (0) + o n (o) 


(3) 


where Q (epl ) , Q p (int), Q p (poly) and Q n (sub) are the excess charges in the 
epitaxial layer, the interfacial layer, the polysilicon bulk, and the 
substrate, respectively with corresponding time constants x p (epi), x p (int), 
tp(poly) and t n ( sub) , and J p (W E ) is the recombination current density at the 
ohmic contact, x = W^. In (1) - (3) we have neglected the recombination in 
the p/n junction space-charge region which can be evaluated separately [13], 
if necessary. 

Having described the cross-sectional microstructure and composition of 
our devices. Figs. 2 and 3, and defined the current components, (1) - (3), we 
will in the next Section analyze the data in order to: (i) separate the 
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measured current into components and (ii) identify the parameters that 
determine the minority carrier transport in the interfacial layer and in the 
polysilicon film. 


IV. EXPERIMENTAL RESULTS AND DISCUSSION 

A. Summary of Results 

Figure 5 shows the measured current-voltage (I-V) characteristics for 
some representative devices from Table I. The hole saturation current 
densities Jpg for all devices used in this study are listed in column IV of 
Table I. The values for Jpg were obtained by determining the electron 
saturation current J n g < 2 x lO" 1 ^ A/cm^ from the doping density and from the 
measured diffusion length in the p-type substrate (L^ > 9 pm) obtained by the 
X-ray method [14] and then subtracting it from the measured saturation current 
Jg. I-V characteristics were measured both on devices with and without the n + 
contact rings; the values for J p g were identical, but the devices without the 
n + rings had a larger series resistance (see discussion below). The I-V 
curves exhibited an ideal J<=exp(qV/kT) dependence for at least 3-4 decades of 
current. Thus no corrections for the non-ideal space-charge-region and 
surface recombination currents were required. The results reported here were 
measured on devices with an area of 62 x 23 mil^ in which the perimeter 
effects are negligible. Each J p g value in Table I is the average value over 
at least 20 devices. 

The following important observations can be made from Fig. 5: 

(a) The current density for the in-situ doped polysilicon devices is 
about 3 orders of magnitude smaller in comparison with that of the 
metal contact. This demonstrates the sensitivity of our test 
structure to the interface conditions and the usefulness of the 
polysilicon contact in suppressing the current. 
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(b) The heat treatment after the deposition of the in-situ doped 
polysilicon film causes almost no changes in the measured 
current. This is evident from the comparison of the device 1A 
without any heat treatment with the device IB heat treated at 
800° C, 64 hrs. 

(c) The undoped polysilicon film, device IE, behaves as a surface with 
large recombination velocity with no improvement in comparison with 
the reference metal contact. 

(d) The bilayer device II without heat treatment is identical to the 
undoped device IE without heat treatment. 

(e) The bilayer device 1J after 800°C, 64 hrs heat treatment, exhibits 
about the same current density as the in-situ doped device IB with 
the same heat treatment. 

(f) Devices, without the n + contact ring, with undoped film (IE) and 
bilayer film without heat treatment (II) have a large series 
resistance R$, which is demonstrated by the deviation from the 

J = exp(qV/kT) dependence at V > 0.3 volts. The series resistance 
R^ is greatly reduced for these devices if the measurement is taken 
on the device with the n + contact ring, as illustrated in Fig. 5 
for the device IE. 

The above observations (a) - (f) and the results in Table I lead us to 
the following conclusions: 

(1) Large currents measured in all devices with the undoped polysilicon 
layer contacting the monosilicon layer preclude the possibility of 
an effective thin interfacial oxide layer [2,7] between the 
polysilicon and the monosilicon in our devices. It was proposed by 
de Graaf and de Groot [2] that a 20 - 30 A thick interfacial oxide 
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layer represents a tunneling barrier for the minority holes, which 
suppresses the minority-carrier current in the polysilicon. If 
such a layer was present in our devices, the hole current would be 
limited by the tunneling mechanism [2,5,8] through this layer in 
both the doped and undoped samples. In this work, we did not 
intentionally grow an interfacial oxide prior to the polysilicon 
deposition that may produce such a tunneling barrier [2,7]. 

(2) The polysilicon-monosilicon interface is very sensitive to the 

arsenic concentration, rather than to the electron concentration. 

This conclusion is based on the observed insensitivity, discussed 

in (b) above, of the hole current to the heat treatments of the in- 

situ doped devices. The electron carrier concentration in the CVD 

deposited polysilicon is a function of the postdeposition heat 

treatment [11,12] and increases with the heat treatment 

temperature. Thus, large differences in the carrier concentrations 

can be expected before and after the heat treatments [11,12]. On 

the other hand, the arsenic dopant concentration is almost 

independent of the heat treatment in the in-si tu doped films [9], 

The correlation between the arsenic concentration at the 

interface on the one hand, and the current density and on the 

other hand, is further demonstrated in Fig. 6 for the devices with 

bilayer polysilicon films. The comparison of device II with no 

heat treatment with two devices after heat treatment (devices 

1J,1L) shows a drastic reduction in the current density. The 

arsenic concentration at the interface before the heat treatment is 

18 -3 

the residual <10 cm , while after the heat treatment the 

/V 

20 -3 

concentration is ~10 cm" (see Fig. 2). Device IK heat treated 
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at 900°C for 15 min still gives J p q equal to that of the device 

II. Apparently, the arsenic does not reach the polysilicon- 

monosilicon interface after the 900°C, 15 min thermal treatment. 

This is consistent with the recently published values for arsenic 

diffusivity in the grain boundaries [9]. From Fig. 3 in [9] we can 

deduce that the product (Dt)*/ 2 , w nich is a measure of arsenic 

penetration through the grain boundaries, has to be much larger 

than about 500 A for arsenic to penetrate through the ~ 1500 A of 

the originally undoped layer of the bi layer film. For 900°C, 15 

min heat treatment, the value of (Dt) 1 / 2 is only about 550 A, but 

is about 2000 A for the devices 1J and 1L, where arsenic did reach 

the interface, as Fig. 2(b) illustrates. The device 2BA, for which 

(Dt) 1 / 2 = 280 A, behaves similarly to device IK (see Table I). The 

R$ is large for the non heat treated case and decreases with 

increasing (Dt)^ 2 , thus directly reflecting the penetration of the 

highly doped and conductive arsenic front through the grain • 

boundries towards the interface [9]. A large in all undoped 
18 3 

devices (< 10 As/cm ) and in the bi layer device II suggests 
that the residual arsenic in the undoped films is present mostly in 
the nonactivating lattice sites, or more probably, that most of the 
free electrons are trapped at the grain boundaries leaving the 
grains essentially depleted [15,16], which results in a large 
resistivity of the polysilicon film. The reduction in R s after the 
heat treatment can also be attributed to the lowering of the grain 
boundary barriers in the polysilicon film because of the increased 
doping in the vicinity of the grain boundaries [16,17]. 
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B. Separation of Current Components 


In the previous discussion we described the electric behavior of devices 
with various polysilicon films and their correlation with the microstructural 
composition. Below, we analyze the data and attempt to identify the possible 
dominant physical mechanisms governing the minority carrier transport in the 
polysilicon film. In order to do that, we have to first separate the measured 
currents into components defined in (1) and (2) and, thus, extract the desired 
values of the recombination currents in the polysilicon. Then we correlate 
the phenomenological model with experimental results for polysilicon films 
with different microstructures. 

The five current components in (1) - (2) can be divided into two groups, 
based on the criteria of the doping density in the particular region. Regions 
with doping densities smaller than about 10*® cm - ® are expected to exhibit 
activation energy corresponding to the energy gap of intrinsic Si, while 
regions with doping much larger than 10*® cm - ®, will have smaller activation 
energies (bandgaps) because of the bandgap narrowing aEq [18]. Considering 
the above criteria, we can write 


Qp(epi) Q n (sub) -E GI /kT 
Tp(epi) + x n (epi ) a e 


Qp(int) 


Qp(poly) 

Tp(poly) 


-(E 

+ J p (W E ) a e 


SI 


- AE G )/kT 


(4) 

(5) 


where Eqj is the intrinsic Si bandgap, aE g is the average bandgap narrowing 
for the entire heavily-doped region, Q po and Q no are the preexponential 
constants in Q - Q Q exp(qV/kT) , and J p0 (W E ) is the saturation current 
density. The activation energy in (5) is thus Eg = E GI - AE y < E G j . As 
discussed earlier, the precise value of the Q n0 (sub)/x n (sub) component was 
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already determined independently, thus, only the hole current will be 
discussed henceforth. 

Figure 7 shows the activation plots [18] for representative devices, with 

undoped, bilayer and in-situ doped polysilicon films. The activation curve 

for the reference device with metal contact is also shown. 

To construct the activation plot, we have to properly account for all 

temperature dependent parameters in (4) and (5). This is done by plotting 

[18] C(Jp 0 /T m )[exp(n c )/F 1/2 (ric)] 3 exp[-(E GI ( 0 ) - AE G )/kT] = exp(-E G /kT) 

versus 1/T . Here, E G j( 0) = 1.206 eV is the linearly extrapolated energy gap 

of intrinsic Si at 0 K, C is a temperature independent constant that involves 

device parameters (such as geometry, doping, minority-carrier mobility, and 

minority-carrier lifetime), n c = (Ep N - E c )/kT (see Fig. 4(b)), and F^ 2 (r! c ) 

is the Fermi -Dirac integral of the order of 1/2. In the above expression, m=4 

for a transparent region bounded by a surface with very large surface 

recombination velocity, e.g. metal contact or undoped polysilicon contact, 

2 4 

where Jp G “[n^ (kT/q)^] « T , provided that the minority-carrier mobility pp 

is temperature independent [18]. On the other hand, if the metal contact is 

replaced by a contact with low surface recombination velocity, then 
2 3 

OpO «(n.W/Xp) « T [19], i.e, m=3, provided that the minority-carrier 
lifetime is temperature independent. 

In Fig. 7, we show the activation plot for the in-situ doped device IB 
and the bilayer device 1J for which Jp Q (polysilicon) << J pG (metal); thus, 
m=3 . This yields E Q = 1.206 t 0.03 eV (i.e. aE g = 0 ± 30 meV) . The 
assumption of Xp / f(T) is very reasonable in the narrow temperature range, 
300 K < T < 370 K used in Fig. 7. The slope of the activation plot is 
relatively insensitive to small variations in m. For example, the activation 
energy for the device 1J for m=3.5 is larger by only about 10 meV compared to 
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that for m=3. Thus, only a negligible error is caused if x p varies slightly 
with T. 

In the case of the devices with metal contact and undoped polysilicon 
contact films, setting m=4 neglects the temperature variation of the minority- 
carrier mobility p p . The result for m=4 is Eg = 1.146 ± 0.020 eV, i.e. AEg 
* 60 ± 20 meV. Assuming now that the temperature dependence of the minority- 
carrier hole mobility in the low-doped (~10^cm“^) epitaxial layer is the same 
as the temperature dependence of the majority-carrier hole mobility, which 
follows Up a T - *’^ for 300 K < T < 400 K [20], then a correction to the above 
result is required. This is done in Fig. 7 using t^T - ^*^ = T^'^ instead of 
T 4 , yielding Eg = 1.206 ± 0.025 eV, I.e. aEq = 0 ± 25 meV. 

The activation energies for the other devices from Table I, not displayed 
in Fig. 7, are equal to the values of the representative devices discussed 
above. 

In contrast. Fig. 7 also illustrates the activation curve corresponding 

to the hole recombination in the n + monosilicon layer obtained from the n + /p 

monosilicon junction with uniformly doped transparent 0.15 pm n + layer with 

19 -3 

electron concentration N = 2 x 10 cm . The activation energy for this 

n + -monosil icon layer (using m=4) is smaller by about 0.135 eV, Eg * 1.071 ± 

0.025 eV, in comparison with the value for the epitaxial n-monosilicon 

obtained from the device with the metal contact. This difference can be 

interpreted as bandgap narrowing in the n + -monosil icon region of 

AE„ * 135 meV [18]. 
b 

The activation energy of the reference device with the metal contact is 
described exactly by (4) (polysilicon is absent here) and the current is due 
completely to the low-doped epitaxial layer (J po >> J no ) where AEg = 0. The 
fact that the result obtained from the activation plot of this device gives 
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E g = EqjCO) = 1.206 eV within the experimental accuracy of about + 0.026 eV 
establishes a credibility for the activation plots and also provides a value 
for comparison with the devices with polysilicon contacts. 

The hole saturation current density for the devices with polysilicon 
contact in which J po (polysilicon) « J po (metal), assuming L p » W ep -j and S p << 
Dp/W ep i, can be written as [19,21] 


qn. 


PO “ N 


DO 


epi 

-Tp(epi') 


V 


(6) 


Here S p is the effective surface recombination velocity at x = W epi , resulting 
from the recombination losses in the polysilicon. The observed activation 
energy and the fact that J pQ (polysilicon) « J p0 (metal) imply S p « 

W -/t (epi). The hole distribution corresponding to this case is 
schematically depicted in Fig. 4(c). From (6) we then find T p = 0.6 psec and 
L p = 26 pm. 

To corroborate the fact that S p is negligibly small in our devices and to 
avoid the uncertainties in the interpretation of the activation plots, we have 
fabricated the test structures by depositing an in-situ doped polysilicon 
films on epitaxial layers with different values of W epi : 0.8 pm (p = 0.6 Qcm) , 
8 pm (p = 1.2 Qcm) , and 12.8 pm (p = 0.6 Qcm) . The deposition and the heat 
treatment conditions were identical for all three devices. The heat treatment 
was done at 900°C for 15 min, and corresponds to that of the device 10 from 
Table I. 

Using (6), we can calculate the values of S p + W epi /t p (epi) from the 
measured J pQ and from the known doping densities. Then a plot of S p + 
w epi / T p^ ep i ) versus w epi is ex P e ^ed to yield a straight line, with a slope 
proportional to T p (epi) = 0.6 psec, assuming that t p in the 8 pm and 12.8 pm 
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epitaxial layers is about equal to that in the 0.8 pm layer deduced above. 
Since Sp is identical for all devices, an extrapolation to 0 yields 

S p . This plot is shown in Fig. 8, yielding S p ~ lb cm/sec « W epi /T p (epi) for 
all three devices. 

The exact value of Sp is of importance, since it represents the desired 
recombination losses in the polysilicon, including the losses in the 
interfacial layer. The approximate value of S p ~ 15 cm/sec deduced above 
indicates that the recombination losses in the polysilicon film are only about 
10fc of the measured hole current. The polysilicon-monosil icon interface thus 
acts as a very efficient barrier for the minority carriers. 

Hence 


Tplpoly) - + 


<W int > 

•tp(i nt) 



A/cm^ 


(7) 


where we have neglected Jp 0 (Wg;) because it was found experimentally [3] that 
Lp(poly) < 500 A, which is much less than the thickness of the polysilicon 
films used in this work. 


C. Discussion 

The TEM study revealed that the grains in the in-situ doped polysilicon 
films after heat treatment at 900°C - 1000°C are essentially columnar and the 
grain diameter is at least 400 A and as large as about 1000 A [22]. The 
effective hole diffusion length in the heavily-doped polysilicon is Lp < 500 A 
[3], i.e., Lp < grain diameter. Such polysilicon film thus has only a few or 
no grain boundaries in the direction of the current flow and can be considered 
to be almost a single crystal n + layer separated from the monocrystalline 
epitaxial layer by only grain boundaries (or an interfacial layer) at the 
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metallurgical interface. Thus, if the very thin (~200 A) interfacial layer 
(Figs. 3, 4) is transparent for the minority holes, then the diffusion hole 

IQ O 

current in the bulk of the polysilicon film with N ~ 10 cm is expected to 
be about equal to the hole current of the single crystal n + /p junction of 
Fig. 7 with N = 2 x 10^ cm ^ which is Jp 0 = 4 x 10 ^ A/cm^ at 25°C. This 
value is about one decade larger than the measured J p0 for the devices with 
the in-si tu doped polysilicon films and at least a factor of 40 larger than 
the value estimated in (7) for the total recombination losses in the 
polysilicon film. To explain the discrepancy of about two decades in the 
current, the effective transport parameters in the polysilicon would have to 
be markedly different from their monosilicon values. 

Modeling of the transport in the polysilicon is complicated by the grain 
boundaries. As a first approximation, we will treat the polysilicon 
analogously to the transport in the single crystal Si, as justified above, but 
with an effective value for the parameters, such as bandgap, mobility, and 
lifetime, to account for the influence of the grain boundaries. The effective 
parameters are average parameters for the polysilicon film corresponding to 
the average grain size. This approach was used before [3,4]. We can then 
write 

Qp(poly) q"i o (poly) y 

TpfpWr ‘ ypoTyy CpOpVTy) ex P<ff> • < 8 > 

where N e (poly) * 10 18 cm* 3 is the effective doping [18,23] in the heavily 
doped polysilicon. In (8), we have used the quasi-equilibrium assumption of 
an almost flat hole quasi -Fermi level in the epitaxial layer up to the 
metallurgical interface, Fig. 4(b). This is justified by the fact that 
J po (polysilicon) « J po (metal) and by noting that the recombination losses in 


3-17 



the very thin (< 200 A) diffused region in the epitaxial layer can be 
neglected [19.24J. Note here that, to remove the discrepancy, Dp(poly) would 
have to be at least 40 times smaller than the single crystal value. This is 
not a likely possibility in the large grains with Lp ~ grain diameter. To 
explain the results by lower value of aEq in the polysilicon requires a 
difference of about 100 meV between the single crystal and polysilicon 
values. However, the measured activation energies in our devices and in the 
bipolar transistors [3] do not support this assumption. 

For the bi layer devices, we have a more complicated situation. The 
interfacial region is heavily doped, followed by a 1500 A thick undoped layer 
and again a heavily doped 1000 A-wide in-situ doped layer (Fig. 3). The 
heavily doped grain boundaries (N ~ 10 20 cm -3 ) within the lightly doped grains 
create a high-low junction with low effective recombination velocity [25], 
thus only slightly affecting the current flow. 

If the minority holes were to diffuse across the about 200-300 A thick 
interfacial layer into the low-doped region, the expected current from (8) 
would be much larger than for the in-situ doped devices, because N g (poly) << 

1 Q O 

10 10 cm in the low-doped grains. The measured hole currents for the in-situ 
doped and bilayer devices are, however, almost equal. The exceptions are the 
bi layer devices where the heat treatment was not sufficient to obtain high 
arsenic concentration at the interface. The differences in the values of Jp Q 
of less than about a factor of two between the bi layer devices 1J, 1L and the 
in-situ doped devices, can be attributed to small variations in , -Up(epi), 
and mainly N^. The average doping density in the epitaxial layers was 
measured by the C-V technique and varied from about 9 x 10 cm“ J to about 
2 x 1016 cm" 3 . This fact is also reflected in the spread of the Jp 0 values 
for the devices with the metal contact which was in the range from about 
4 x 10“10 A/cm 2 to about 6 x lO"^ A/cm 2 . 
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The preliminary conclusions (1) and (2) in Section I V-A and the 
discussion in this section lead to the following: 

(i) The minority carrier transport at the polysilicon-monosilicon 
interface with no intentional insulating layer at the interface is 
determined by the carrier transport within the polysilicon, rather 
than by a tunneling through an insulating layer. 

(ii) The similarity of the electrical properties of the devices with the 
in-situ doped and the bilayer polysilicon films suggests a 
correlation between the electrical properties and the 
microstructure of the very thin interfacial layer within about 

20U - 300 A of the metallurgical interface in the polysilicon 
film. It also suggests that most of the minority carriers crossing 
the x = W epi plane (Fig. 4) recombine within 200 - 300 A of the 
metallurgical interface. 

(iii) In order to explain J p (polysilicon contact) « J p (n + monosilicon), 
the heavily-doped interfacial layer must have properties 
significantly different from the properties of the monosilicon. 

(iv) The dominant mechanism responsible for low J p (polysi 1 icon contact) 
is much more sensitive to the arsenic dopant concentration at the 
interface than to the carrier electron concentration at the 
interface. 

The polysilicon-monosilicon interface can be regarded as a region with a 
high degree of disorder. This is due largely to the very sudden transition 
from the orderly epitaxial single crystal layer to the CVD-deposited 
polysilicon with the randomly oriented (without prefered orientation) grains 
of a few hundred Angstroms in size. It was recently observed by means of 
cross-sectional TEM that the arsenic doped polysilicon-monosilicon interface 
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consists in some cases of a crystalline phase with about 40 A average grain 

size, compared to ~ 400 A in the bulk of the polysilicon film, while in some 

other cases this interface is amorphous [22,26]. It is well known that the 

disordered regions of semiconductors are associated with low carrier 

mobilities (and diffusivities) . As an example, the drift majority-carrier 

-2 2 

mobilities in amorphous Si are a (2-5) x 10 cm /Vs and * 5 x 
10"^ cm^/VS [27]. Even smaller values are expected for minority carriers. 

The disorder contributes to the random component in the atomic potential 
which results in band tails [28]. Trapping by band tails was identified as 
one of the possible mechanisms causing very low mobility in amorphous Si 
[27]. Band tails can be responsible for low minority-carrier mobility by a 
process of trapping also in heavily doped single crystal Si, as was recently 
suggested [29]. 

The microstructure at the interface and the electrical measurements lead 
us to a proposal for a possible mechanism of current suppression by the 
polysilicon contact: very low minority carrier mobility (and diffusivity) in 
the very highly doped and disordered interfacial region. Note that both 
amorphous structure and very small grains at the interface predict low 
mobility. 

Because of the large sensitivity of our test structure to the interface, 
we can make a better estimate of the minority-carrier mobility in the 
polysilicon than made before in [3], In order to do that we have to have a 
reasonable model for the recombination current Q (int)/x (int) in the 

H r 

interfacial layer. We assume again that, as a first approximation, this 
current can be described as a diffusion current 


Qp(int) qn* 
'Cp('int) = N e (int) 


Dp(int) 

tp(inty 


< 10 


-13 


A/cm^ 


(9) 
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The quantity on the right-hand side of (9) results from (7). Approximating 
the parameters in (9) by a single crystal Si values, we have 
N g ~ 10*^ cm~^ [18], and Tp(int) ~ x^ = 2 x lO'^sec [30], where x ^ is the 
Auger lifetime. Equation (9) then implies the upper limits for 


Dp(int) ~ 5 x 10~ 3 cm^/s 


(10a) 


and 


Lp(int) ~ 100 A . (10b) 

Thus, Op(int) from (9) and (10a) is about 3 orders of magnitude smaller than 
the single crystal value. 

The rough approximations of (10a), (10b) do not take into account the 
consequences of the detailed microstructure of the interface, but they do, 
nevertheless, lend support to the interpretation based on the possibility of 
large gradient in the minority-carrier diffusivity at the polysilicon- 
monosilicon interface. This mechanism will greatly suppress the current. The 
minority carriers will then recombine very close to the metallurgical 
polysilicon-monosilicon junction without the possibility to diffuse deep into 
the polysilicon. 

An alternative explanation can be offered by an interface in a 
heterojunction form with a very large barrier for holes. The heterojunction 
can be created because of a position-dependent chemical composition giving 
rise to spatial variations of the bandgap, electron affinity and the band 
edges of the conduction and valence bands [31]. As an example, if a thin 


3-21 


(~ 20 A) amorphous layer exists at the interface [22,26], a discontinuity of 
the bandgap and electron affinity may result. To further analyse such a 
structure would, however, require much more detailed information about the 
chemical and electronic structure of the interface than presently available. 

Even though the heterojunction model or the tunneling model [2] cannot be 
excluded as a possibility, there are however two findings that point 
particularly in the direction of the low mobility model. First, the finding 
of the correlation between the recombination current and the dopant arsenic 
density agrees with studies which show that at the microscopic level the 
atomic potential has a random component because of the random distribution of 
impurities whether impurities are ionized or not [32]. Second, the 
observation by Ning and Isaac [3] of the dependence of the base current in the 
bipolar transistor on the polysilicon thickness. A noticeable increase in the 
base current was observed for transistors with the polysilicon film thickness 
of 300 A or less. This observation agrees with our model of carrier 
recombination within the interfacial region, such as illustrated in Fig. 3 for 
the bi layer film structure. 

The undoped polysilicon contacts in devices IE - 1G behave as interfaces 
with very large recombination velocity with current densities about equal to 
the current densities for the reference metal contacts (see Table I). This is 
expected because the potential profile at the interface between the epitaxial 
layer with N = N Q0 = 10^ cm”^ and the polysilicon with N < 10 ^ cm“^ creates 
an electric field in the direction of the hole flow, thus aiding holes across 
the interface. The recombination current inside of the undoped contacts is 
expected to be large because of low electron density and small grain size 
(~ 100 A), as (9) implies. The above arguments apply also to the bi layer 
devices II and IK. 
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It is very convenient to characterize the polysilicon-monosilicon 
interface by the effective surface recombination velocity S p at the high-low 
junction in the epitaxial layer at x = W epi (Fig. 4), which is defined as [33J 


S = J (W .)/qAP(W .) 
p p x epv M v epr 


(ID 


where AP(W epi ) is the excess hole density and J p (W epi ) is the current density 
at x = W epi . The value of S p cannot be deduced exactly without a precise 
knowledge of T p (epi), since the bulk recombination in the epitaxial layer is 
not negligible in our devices, but is in fact dominant. An approximate value 
of S p ~ 15 cm/sec was deduced before. For many purposes, however, an upper 
limit of S p is sufficient. This upper limit can be obtained exactly from (5) 
by neglecting the bulk losses (i.e. by setting W epi /-c p (epi ) =0), or from 


pmax 


V%i 

T3757737 


( 12 ) 


where J p (Q) = [(qn^D J/INqdW . )] exp(qV/kT) is the current for the reference 
device with the metal contact. The upper limit values for S p calculated from 
(12) are listed in Table I. 

From the upper limit of S p given by (12) we can also estimate the upper 
limit of the effective surface recombination velocity at the n + polysil icon-n + 
monosilicon metallurgical interface. Assuming that the thin n + outdif fused 
monosilicon layer is essentially transparent we obtain from J p (W epi ) - 
J p (n + poly-n + ) 


S p (n + poly-n + ) » S p 


N e ( int) Op(int) 
N 0D (epi) = 7 Tp(i'ntJ < 10 


cm/ sec 


(13) 
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where N 0 ( i nt)/N DD (epi ) = 10 18 cm“ 3 /10 16 cm" 3 = 100 [18]. Using the 
approximate value of S p ~ 15 cm/sec in (13) yields S p (n + poly-n + ) = 

1500 cm/sec. 

The n + polysilicon contact on the top of the n + diffused or implanted 
monosilicon emitter in the bipolar transistor can be also characterized by 
(13), since (13) reflects the recombination in the polysilicon independently 
of the doping in the monosilicon. The low value for S p (n + poly-n + ) in 
comparison with S p ~ 10^ cm/sec for the ohmic metal contact explains, 
qualitatively, the improved current gain in bipolar transistors with 
polysilicon contacts to thin emitters. 

IV. SUMMARY 

In this work, we have presented the results of a detailed experimental 
study of a polysil icon-monosilicon contact. The purpose of this study was to 
explain both qualitatively and quantitatively the mechanism of improved 
current gain in bipolar transistors with polysilicon contacts. The study 
correlates the microstructure of the polysilicon-monosilicon junction with the 
electrical properties. The electrical properties were measured using p/n 
junction test structures (Fig. 1) that are much more sensitive to the contact 
properties than are bipolar transistors. 

Correlation of the electrical results with the microstructure yielded the 
conclusion that the minority-carrier transport in the polysilicon-monosilicon 
junction is controlled by a highly disordered layer within ~ 100 A of the 
interface characterized by very low minority-carrier diffusivity. This 
conclusion confirms the model of Ning and Isaac that explained the improved 
current gain by D p (n + polysilicon)/ D p (n + monosilicon) = 0.3. The present 
results indicate, however, that the effective D p in the n + polysilicon has to 
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be substantially smaller than the corresponding value in the n + monosilicon to 
obtain the quantitative agreement with the data: D (n + polysilicon)/D (n + 

r r 

inonosilicon) ~ 5 x 10"^. The effective surface recombination velocity for 
minority holes at the n + polysilicon-n + monosilicon interface was found to be 
much smaller (~ 10^ cm/s) than the kinetic limit value for ohmic contact 
(~ 5 x 10® cm/ sec). 

We did not study the effects of surface treatments prior to polysilicon 
deposition, namely the effects of thin interfacial oxide layer [2,7]. 

The saturation current density in the polysilicon contact was found to be 
less than about 10 A/cm at 25°C. The typical value of the saturation 
collector current density in the bipolar transistor with basewidth of about 
0.1 pm, base sheet resistance of about 8 kQ/square and current gain p = 100 is 
about 10“ u A/cm [3]. This implies that current gains p > 1000 are possible, 
provided that the recombination losses in the n + monosilicon emitter can be 
suppressed, for example by reducing the width of the n + monosilicon region 
[2,34], 

The sensitivity of the test structures used in this study to the 
polysilicon contact and the accuracy of the method used to separate the 
current components was limited mainly by the minority-carrier lifetime in the 
epitaxial layer. More conclusive experiments [21] using devices with varying 
widths and dopings of the epitaxial layer would be desirable to improve the 
accuracy of the analysis. 

A very low effective surface recombination velocity for the minority 
carriers at the polysilicon-monosilicon interface suggests applications of 
polysilicon contacts in devices other than bipolar transistors. An example 
could be a back-surface-field solar cell in which the recombination losses at 
the back high-low monosilicon contact can limit the cell efficiency. 
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TABLE I 


Summary of heat treatments and electrical measurements. Heat treatments for 
poly-silicon contacts were performed in argon. Electrical measurements were 
done at 25°C. 


CONTACT 

DEVICE 

HEAT 

TREATMENT 

^po 

(10"^A/cm^) 

S 

pmax 

(cm/sec) 

Undoped 

IE 

none 

(5.8 ± 0.4) x 10 3 

~ 10 6 

Polysilicon 

IF 

800°C, 64 hrs 

(5.5 ± 0.5) x 10 3 

~ 10 6 

(1500 A) 

1G 

900°C,-5 min 

(5.5 ± 0.5) x 10 3 

~ 10 6 

In-situ 

1A 

none 

7 ± 0.3 

175 

Doped 

IB 

800°C, 64 hrs 

10 ± 1.5 

250 

Polysilicon 

1C 

900°C, 5 min 

6.6 ± 0.5 

175 

(1500 A) 

ID 

900°C, 15 min 

6.5 ± 0.3 

165 


2BE 

1000°C, 15 min 

6 ± 0.5 

165 

Bilayer: 

11 

none 

(4.2 ± 0.2) x 10 3 

~ 10 6 

1500 A undoped + 

1J 

800°C, 64 hrs 

4.9 ± 0.6 

175 

1000 A in-situ 

IK 

900°C, 15 min 

(4.2 ± 0.2) x 10 3 

~ 10 6 

doped 

1L 

850°C, 14 hrs 

5.2 ± 0.6 

180 


2 BA 

750°C, 8 hrs 

(4.0 ± 0.3) x 10 3 

~ 10 6 

Metal 


450°C, 20 min 

(5 ± 1) x 10 3 

~ 10 6 

Reference 
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FIGURE CAPTIONS 


Fig. 1 Cross-sectional view of devices used in this study. 

(a) Device with polysilicon contact film; 

(b) Reference device with metal contact only. The broken lines in (a), 
(b) show the implanted n + -contact ring in the epitaxial layer. 

Devices with and without the contact rings were fabricated on the same 
wafer. 

Fig. 2 SIMS profiles of the arsenic distributions [8] for: 

(a) The in-situ doped polysilicon with no heat treatment; 

(b) The bi layer film before heat treatment (II) and after 750°C, 8 hrs 
(2BA) and 850°C, 64 hrs (1J) heat treatments. The extent of the As 
profile into the monosilicon is not larger than the resolution of the 
SIMS equipment which is about 1 decade per 100 A [9j. 

Fig. 3 Schematic illustration of the arsenic concentration profile in the 

vicinity of the polysilicon-monosilicon metallurgical interface for the 
bi layer film after heat treatment. The shaded areas show the 
preferentially doped regions around the grain boundaries and the arsenic 
diffusion into the monosilicon. The broken lines depict the width of 
the arsenic diffusion from the grain boundary monolayer to the bulk of 
the grains. The screened area within the interfacial layer symbolizes 
the fact that the microstructure of the polysilicon film in this area is 
very different from the microstructure in the rest of the polysilicon 
film. 

Fig. 4 (a) Schematic illustration of a forward biased test structure with 

heavily-doped n + polysilicon contact showing the separation of the 
polysilicon film into two layers: the interfacial layer and the 

polysilicon bulk. 

(b) Band diagram for the forward bias. 
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(c) Profile of the minority hole concentration for the case of low 
effective recombination velocity at x = W e p^ . The broken line depicts 
the approximate decay of the minority holes in the polysilicon. 

Fig. 5 Current-voltage dependencies for several devices with polysilicon 
contact and for the reference device with metal contact. 

Fig. 6 Current-voltage characteristics for devices with bilayer polysilicon 
contact . 

Fig. 7 Activation plots [18J of C( I p0 /T m )Cexp(rj c )/F 1/ 2C T ic^ ^ = exp[-(E GI /(0) - 
AE G )/kT] = exp(-E G /kT) versus 1000/T for devices with undoped, in-situ 
doped, and bilayer polysilicon contacts. Activation plots for the 
device with reference metal contact and for the reference n + /p 
monosilicon junction are also shown. The slopes yield Eg * 1.206 ± 0.03 
eV for the metal reference device and for the polysilicon devices 18, 

1J, and IF. The slope for the n + /p monosilicon device yields Eg = 1.071 
± 0.025 eV. 

Fig. 8 The dependence of S p + W epi /xp(epi) on the thickness of the epitaxial 
layer W . . The extrapolation to W ^0 yields S ~ 15 cm/sec. 
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CHAPTER FOUR 

SOLUTION OF THE CONTINUITY EQUATION 
IN PLANAR SYMMETRY CASES AND ASSESSMENT 
OF PHOTOLUMINESCENCE DECAY 

I. INTRODUCTION 

Photoluminescence decay is one of the methods used in determining the 
minority carrier lifetime in heavily doped silicon [1]. The method consists 
of exciting a silicon sample by a laser pulse and thus generating excess 
electron-hole pairs. The lifetime is derived from the decay rate of the 
excess charge. This decay as a function of time is obtained by monitoring the 
radiation emitted by the few electron-hole pairs that recombine through 
radiative recombination. In principle, this method is more reliable than many 
other methods because it involves no assumptions regarding the band structure 
and the transport parameters in heavily doped silicon. On the other hand, the 
method has the drawback that surface recombination prohibits the decay from 
being purely expoential. Moreover, the effect on the decay rate of the size 
and shape of the laser beam was not clarified previously. 

Recent papers [2,3] have investigated the time dependence of the 
luminescence decay, and in particular the influence of surface recombination 
on the decay rate for special types of excitation. In [2] the continuity 
equation is solved analytically for the total excess charge in the case of a 
finely focused beam, while in [3] the problem is simplified to one dimension 
by considering uniform excitation. In both cases it is found that the effect 
of the surface recombination can be minimized by waiting several lifetimes, 
after the laser is turned off, before measuring the decay rate. 

In this work we extend the analysis of [2,3] for the general case of 
arbitrary excitation, and prove that the beam size, shape and uniformity does 
not affect the decay rate. An upper bound for the error due to surface 
recombination is derived. In agreement with [1,2,3J the effect of surface 
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recombination is found significant only in the initial part of the 
transient. In addition, as a preliminary result, a simplified form is given 
for the continuity equation when planar symmetry exists. 

II. SOLUTION OF THE CONTINUITY EQUATION IN PLANAR SYfUETRY CASES 

Figure 1 illustrates an infinitely extending semiconductor sample 
terminating at the z=0 and z=VJ planes. The laser beam shown is meant as an 
example of excitation. The following analysis is, in general, valid for any 
other type of excitation, for example, thermal generation, and minority- 
carrier injection over a potential barrier. The sample is n-type (extension 
to p-type is straightforward). 

The continuity equation for holes is 


r\ p + 

(x,y,z,t) = -V*j p (x,y,z,t) + G(x,y,z,t) - U(x,y,z,t) 


( 1 ) 


where 3 is the hole flux and G and U are generation and recombination rates 
P 

in excess of their thermal equilibrium values. Integration of Eq. (1) over 
the z = Z'(0 < V < W) plane yields 


|2(z',t) = 


/ / V»j dxdy + / / Gdxdy - / / Udxdy 

-CO —os P -OO — OO —CO -OO 


- - 1 I ( 



3j 


3y 


j^-)dxdy - / 


3j 


21 


3z 


dxdy + 


OO OO OO OO 

+ / / Gdxdy - / / Udxdy 
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( 2 ) 
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Here Q(Z',t) = / / P(x,y,Z‘ ,t)dxdy is the number of holes per unit 

OO -CD 

distance on the z=Z' plane, and J px , J p y, J pz are the components 

of 3 = J x + J y + J z. All of the integrals above are assumed to be 

r px \)y yl. 

finite. 

By Gauss' theorem. 


00 00 
J J ( 



+ ~2^-)dxdy = / (j 


px 


X + j 


py 


y) *ds 


where r,,, is the infinite circle on the z=Z' plane, and ds is the 

infinitesimal vector normal to r m . The integral „/ (J x + J y)»ds is 

r „ Px py J 

the hole flux crossing the infinite circle, and is zero under restrictions on 
G(x,y,z,t) that are nearly always obeyed. Therefore, Eq. (2) becomes: 


|| = - J J ■ dxdy + / / Gdxdy - / / Udxdy 


3j 


9Z 


(3) 


where from now on we omit the (infinite) limits of integration for brevity. 

Expressing J pz in terms of drift and diffusion and assuming that the 
electric field, diffusivity, and mobility (E z , D p , y p ) are functions only of z 
and that the recombination rate is linear (U = (P-P q )/t with x = x(z)), we 
rewrite Eq. (3) as 


30 _ n 3 2 Q . 

9t > 3z 2 


3D 

[— 2 - - 
l 9z 


y p E z ] 


90 

3z 


- D 


il^zl 

3z 


1 % 

+ ~] Q + — + // Gdxdy 


(4) 


Here Qq = / / Pgdxdy is the equilibrium number of holes per unit distance on 

the z=Z‘ plane. If, in Eq. (4), P replaces Q and G replaces //Gdxdy, the 
usual one-dimensional continuity equation results. 
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Equation (4) is a differential equation for the total charge on the z=Z' 
plane. It is linear provided E z , Dp, pp are independent of the injection 
level (independent of P) . Assuming this, we can rewrite Eq. (4) in terms of 
the number (per unit distance) of excess holes in the z=Z' plane, 

Q* * / / APdxdy = Q - Q q : 



Equation (5) results from Eq. (4) because -Op aP Q /dz + p p E z P Q = 0 . 

The recombination velocities S 1 , S 2 of the front (z*0) and back (z=W) 
surface of the sample enter to determine the boundary conditions of Eq. (5): 

Sj_ Q'(t.O) = D p |f'(t,0) , S 2 Q‘(t, W) = -D p - d - Q - , - ^ - * W ) . 

Under such boundary conditions the solution of Eq. (5) depends only on 
the total generation rate on the z=Z' plane (JjGdxdy) and not on G as a 
function of x and y. Knowledge of Q'(z,t) does not, of course, provide 
AP(x,y,z,t). However, it does provide the total minority current in the z 
direction, I z = -e Dp SQ'/9z + e Pp E z Q' , and the total number of excess 
holes, 

JJJ APdv = /Vdz 
0 

Therefore, when one seeks to determine either the total minority current 
or the number of excess holes, Eq. (b) suffices, provided only that the 
transport parameters of the semiconductor do not depend on the injection level 
and are functions only of z. Equation (5) is much more convenient to use in 
the determination of (measurable) current than is the standard continuity 
equation. 
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III. PHOTOLUMINESCE OECAY RESPONSE 
We now apply Eq. (5) to the study of the decay of the excess charge 
induced by a laser pulse in a semi -infinite (W=<») sample having constant 
doping concentration. For low injection minority carriers flow by diffusion 
and for D independent of z Eq. (5) simplifies to 


Ml = 

at 


Ol _ 2' 


az 


+ G t (t)(l-R)a e" az 


( 6 ) 


where G t (t) is the total photon flux of the laser beam, R is the reflection 
coefficient, and a is the absorption coefficient. In (6) the subscript p is 
dropped for simplicity. 

The boundary conditions are 


q'(o,z) = o, q'(t,®) = o, Q*(t,0)s = p M . ' . U.P) (7) 


where S is the surface recombination velocity. 

If G t (t) is the unit impulse function, G t (t) = 6(t), then Laplace 
transformation of Eq. (6) with time yields 


5 - (s +h Q' = -(l-R)a e" az 


az 


( 8 ) 


where s is the Laplace transfonn operator. The boundary conditions for 
Eq. (8) are the Laplace transformation of Eqs. (7). The solution is 


Q'(s,z)= 


ae~ aZ ( l-R) 


^ — — -(r+«)a/D(i-R) - M Jl.. e~ /Ls+(1/ ^ J/Q I (9) 

-Da +S+( 1/ x) C-Da +S+(l/T)J[s+(l/T)-0r^] 


where r = S/D. 
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If Q'j-(t) the total number of excess carriers and Q'^Cs) its Laplace 
transformation, then 


Q' t (t) = / Q'(t,z)dz => Q* (s) = / Q'(s,z)dz 

z 0 Z 0 


1 - R 


-Doc +S+(1/t) 


(i-R)(r+g)gD[/s+(i/Ty - /d r] 

/s+(l/ -cT [-Dg 2 +s+( 1/ x)J[s+(l/-r)-Dr 2 ] 


( 10 ) 


Taking the inverse Laplace transform of Eq. (10), we obtain 


Q' t (t) - e-V'ci e^’erfcCU/t/V)- ± e 


2 2 

L r t/x 


erfc(UVt/x)J ( 11 ) 


where L = /Ur is the hole diffusion length. This result holds independently 
of the beam shape and is the same as Eq. (8) of Ioannou et al [2] derived in 
the case of a finely focused beam. 

If the laser pulse has a flux N(t) that lasts for a finite time T, then 
the response Q" t (t) will be the convolution of N(t) with Q' t (t): 


Qj(t) - J N(x') Q^t-x'Jdx* . (12) 


If N(T) is a square pulse of amplitude N and width T then, Eq. (12) 
yields for t > T 




(l-R)rg 

r-g 


-t/x 

N x e 


C(L 2 a 2 -l) 


-1 -1 


g 


2 2 

La t/x 

e erfc(Lg/t/x) - 


T/ L 2 g 2 (t-T)/x 

-e /T e erfc(La/(t-T)/x) 


+ L((L 2 a 2 -l) _1 - (L 2 r 2 -l) _1 ) 
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e T / T e ( t_T )/ T e r f c (/(t- T )/i;) - e t ^ T erfc(/t 7 T) + (L‘r^-l)~^r~* 

2 2 2 2 

e T/x e L r (t-T )/ T er f c (Lr/(t-T)/x) - e L r t ^ i: erfc(Lr/t/x) ] . ( 13 ) 

If S = r = 0, Eq. (13) reduces to 

Q" t (t) = N x (l-R)[e T/T -l] e -t/ \ t > T 

which is a purely exponential function with time constant x. 

As a more general statement, if S=0, the decay of the excess charge will 
be exponential independently of the shape of the surface at which the sample 
terminates and independently of the presence of internal fields. This can be 
proved by integrating the continuity equation. 


5P 1 

at 




P 
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over the volume of the sample and applying the divergence theorem: 


aq' 

at 



Here Q' is the excess mi nority-carrier charge and the integral above is the 
hole flux at the surface. This vanishes for S= 0 . Therefore, 


aQ'_ 
at " 



and Q' 


Q'(0) e' t/x 


t >0 . 


When S^O, eq. (13) implies that Q" t (t) is not exponential. A measure of 
the deviation of Q" t from the exponential function (Constant exp(-t/x) follows 
from considering the function: 
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-dDnQ" t (t)] dQ^/dt _ 1 db ( t ) /dt 
dt ' ' Q" t = x b(t) 

where b(t) is the negative of the expression in the square bracket of 
Eq. (13). The smaller that [db(t)/dt)/b(t)] is the more accurate is the 
estimation of x by the slope of lnQ" t . 

One can prove that [db(t)/dt]/b(t) is an increasing function of S and 
1/T. As a worst case, we consider infinite recombination velocity and impulse 
excitation. Then Q"^. becomes: 

2 2 

Q" t = N'(l-R) e L “ t/x erfc(La/t7x) e" t/x (14) 

where N‘ is the strength of the pulse. In such a case: 


db(t) ,22., 

dt~ _ (l//it) - [La/t/x e L a WT ] erfc(La/t/x) La/t/* 
~b(t) ,2 2 +/ t 

e L a t/x erfc(La/t/x) 


(15) 


This is an increasing function of La and when La -► ®, then 
[db(t)/dt]/b(t)->l/2t. Therefore, in the worst case, by measuring the slope 
lnQ“ t at time t, the time constant derived will be T/[l+(x/2t)] instead of 
x. Thus x/2t is an upper bound for the percentage error in determining x. 
This bound is much lower than the one, 2.47 x/t, considered by Tyagi et al 
[2]. Consequently, the maximum error of Dziewior and Schmid [1] due to 
surface recombination is approximately x/[2(5x)] = lOit. 

If an accuracy of 20% in the recombination lifetime is required, then 
0.8 = l/[l+(x/2t)J which implies t=2x. Now the question arises: how can one 

know when t is equal or greater than 2x? If La/t/x >2 Eq. (14) can be 
simpl i fied: 
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(16) 


-t/T 

Q\(t) - Q" t (0) — * . 

/rcLa/t/x 

This implies that Q" t (2'c)=0.05 Q" t (0) 1/La. If L is not known but one can 
estimate its range, then a can be chosen so that La<5, implying Q" t (2x)>0.01 
Q" t (0). Therefore, by waiting until the response drops to 1% of its initial 
value, t will be greater than 2x and the error less than 20% 


IV. SUMMARY 

The excess charge stored in a sample after laser pulse excitation depends 
only on the flux of photons of the laser beam. The decay rate of the excess 
charge is higher than the reciprocal of the bulk lifetime by, at most, a 
factor x/2t. This is a result of surface recombination, the associated error 
of which can be suppressed below 20% by waiting until the response drops to 1% 
of its initial value when La<5. Thus, in so far as the effects of surface 
recombination and of uniformity of the incident beam are concerned, 
photoluminescence decay is a reliable method for determining bulk lifetimes. 
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FIGURE CAPTIONS 


Figure 1 


Schematic illustration of an infinite semiconductor sample bounded 
by planes z=0, z=W, and excited by a laser beam. 
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CHAPTER FIVE 

TRAP CONTROLLED MINORITY-CARRIER MOBILITY 
IN HEAVILY DOPEO SILICON 

I. INTRODUCTION 

Recently [1], we presented experimental evidence for low diffusivity and 
mobility of minority holes in highly arsenic doped (concentrations 
~ 10 20 cm -2 ) n-type Si. This evidence suggests that the minority-carrier 
diffusivity (and mobility) is about an order of magnitude smaller than the 
majority-carrier diffusivity in comparably doped p-type Si. As an 
explanation, we suggested a simple transport model that emphasized trapping by 
localized tail states in the minority band. The commonly used assumption of 
an equal minority and majority -carrier mobilities and the neglect of the tail- 
states effects on the minority-carrier mobility was also questioned in a 
review article by Abrams et al . [2]. Here we report qualitative evidence that 
minority-carrier diffusion in n + Si may be trap limited. 

Figure 1 illustrates the band structure, including the tail states, for 
n + Si. As was already discussed in [lj, for low-level injection, the minority 
holes occupy energy levels near the mobility edge E v . In the model proposed, 
the holes from the extended states in the valence band can be captured by the 
localized tail -states for some mean time t trapp .j ng , and then released back 
into the valence band. This process, shown by arrows in the Fig. 1, will 
decrease the hole mobility. If the scattering rate of holes inside the band 
is comparable with the hole trapping rate, Sah and lindholm [3] derive that 

fc p = et coll /m p * 1/t coll = 1/l scatt + ^trapping * U) 

Here t co ^-| is the characteristic collision free time dependent on the mean of 
the scattering mean free time, t $cdtt , and the mean trapping time, t trapp -j ng . 
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is the 


of a minority hole trapped at a bound state energy level Ey, and m* 
effective mass for holes. We investigate the temperature dependence of by 
estimating t trspp1 „ 9 [3]: 

'/‘trapping “ e p * A1 " «>> ' E V>' kT 3 (2) 

where A is a temperature independent constant. In (2) we have assumed, for 
simplicity, trapping at one shallow level (close to Ey) in the tail band or at 
an acceptor level in the n + region. Such an acceptor level can come from the 
dopant acceptor (boron in our case) in the p-type substrate of a diffused or 
ion implanted n + /p diode. For a parabolic valence band and no degeneracy, 
m = 2 [3]. However, for band-tail distorted band, the value of m is not known 
and must be determined experimentally, as is attempted below. At low and high 
temperatures, respectively, (2) reduces to 

1/1: coll low T ^scatt ’ ^coll high T ~ ^trapping * ^ 

Combining (1) and (3), we obtain the trap controlled minority-carrier mobility 

^(trapping) “ t(e/mJ)/(AT m )j exp(E ft /kT) (4) 

where E^ = Ey - Ey > 0 is the activation energy of the bound-state level Ey. 

II. EXPERIMENT AND ANALYSIS 

Our experiments provide qualitative support for the hole-trapping 
model. The experiments involved a measurement of the temperature dependence 
of the photocurrent response of both n + /p and p + /n photodiodes. The 
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wavelengths of the incident light used create electron-hole pairs almost 
entirely in the heavily doped n + or p + region. Transport of photogenerated 
minority carriers to the collecting p/n junction and the consequent short- 
circuit current I^q depends on the minority-carrier mobility p and diffusivity 
D [4]. Hence I^q and its temperature dependence supplied a vehicle for 
studying p nd D. 

The devices for which we report findings here were fabricated by arsenic 
implantation into 5 Q-cm p-type Si substrates followed oy 1200°C anneal for 
30 minutes. The resulting n + layer was 1.5 pm deep and the sheet resistance 
was 10 £2/square. The As concentration was ~ lO^cnT^ and nearly independent 
of position x over about 0.6 pm below the surface (x =0); for x > 0.6 pm it 
slowly decreases becoming about 3 x lO^cm"'* at 1 pm. The surface was 
unpassivated, covered only by a thin (~ 10 A) native oxide layer. Other 
devices with both unpassivated and passivated surfaces were also 
studied. The results are similar to those to be presented here. The metal 
contacts covered less than 10% of the area of the front surface. 


The I$c(T) dependence was measured in the wavelength range X = 0.38 pm to 
0.4 pm for which the electron-hole generation rate follows G(x) = G(0)exp(-ax) 
[4], where a > 10^ cm -1 [5]. Hence contributions of G(x) to I sc originating 
in the p/n junction space-charge region and in the p-type substrate are 
negligible. Thus [4] 


AqF ( 1 


-Jp 2 - sinh (£-) + cosh (£-) 
P P P 


Here A is the device area, F is the illumination flux density, R" is the 
reflection coefficient, Sp is the effective hole surface recombination 
velocity, W = 1.3 pin is the thickness of the quasi-neutral n + layer, = 
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(Op-cp) 1 / 2 is the hole diffusion length, x p is the hole lifetime, and 0 p is the 
hole diffusivity. 

The analysis of (5) is complicated if one wishes to obtain the magnitude 
of the diffusivity Dp or mobility p.p = Dp(kT/q) from the measured 1^. In 
this note we are interested, however, only in a qualitative trends and try to 
present a simple picture consistent with the data. For this purpose we 
consider now the temperature dependence of 1^. The temperature dependences 
of a and R [6] and -c = x^ (Auger lifetime) [7] are very small. This leaves 
Sp and Dp as the only temperature dependent parameters in (5). We neglect the 
temperature dependence of S p , for a moment, later justifying this assumption 
based on experimental results and on theoretical predictions. 

Recent measurements indicate that L p = 1 pm in the n + silicon doped at 
about 10 20 cm -2 C8]. Using W/L p * 1 we expand the hyperbolic functions in (5) 
and set sinh x - cosh x: 


*SC = ^ f aD p + V (6) 

where B = AqF(l-R)/W is a temperature independent constant. Consider (6) for 

two special cases. First, if S p >> ccDp, =* 8/a. In this case •sc * f (T>- 

This contradicts our experimental result that Icr increases with T (Fig. 2). 

5 -1 

Thus the second case prevails: aO >> S (note that a > 10 cm ). 

P P 

To obtain the simple picture, we try the assumption, m=0, and find from 
(4) and (6) that 


I sc /T » C exp(E A /kT) 


( 7 ) 
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where C is also nearly temperature independent. From (7), log (I$q/T) vs. 1/T 
is predicted to yield the activation energy E A = E-j- - Ey. 

On the other hand, if we make the more common assumption of equal 
majority-and minority-carrier mobilities (and diffusivities) , we would then 
observe 


>SC /T • C ' (8) 

where in deriving (8) we have used the consistent assumption of near 
temperature independence of the majority-carrier mobility (hole mobility = 

45 cm 2 /Vs at T = 20 K; 42 cm 2 /Vs at T = 300 K [9]). 

The experimental data to follow will decide in favor of (7) rather than 
(8), implying trap-controlled minority-carrier mobility as opposed to equal 
mi nority-and -majority-carrier mobilities. 

Figure 2 shows the activation plot of I$q/T versus 1/T for two values 
of F. At low temperature, consistent with our simple model summarized in (7) 
the inferred mobility depends weakly on T because (in our model) the bound 
states in the minority band tail have a negligible emission rate and t scatt is 
a weak function of T. Mobility drops at T > 70 K (where t $catt ~ trapping) 
in agreement with (1). For T > 160 K, p levels off as predicted by (2) for 
shallow traps. From Fig. 2 we can distinguish two temperature ranges E^:= 

10 meV, 70 K < T < 160 K, and E A « 2 meV, T < 70 K. 

These activation energies do not have a simple interpretation because of 
the many levels in the tail. Nevertheless, the value = 10 meV and the 
decrease of pp with increasing T are consistent with (2) and (4) of our model 
for shallow traps with (Ej - Ey) < kT. Moreover, E^ decreases from about 
9 meV to about 8.5 meV as F is increased by a factor of five. Increasing F 
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moves Epp closer to Ey (Fig. 1) which increases the occupation of the tail 
states closer to Ey and E^= (E-p - Ey)/kT is expected to decrease, as 
observed. In the temperature range from about 70 K to about 160 K, where the 
mobility activation behavior is observed, the variation of the energy gap with 
T is very small [10], which justifies our assumption of nearly temperature 
independent a and R [6]. 

The effect of the temperature variation of the effective surface 
recombination velocity Sp(T) on I<jq(T), which was neglected, is now made 
credible by the following considerations. First, if S p (T) is important, E A is 
expected to be about Eq/ 2 = 0.55 eV in a Shockely-Read-Hall model even though 
the surface has distributed states in the energy gap [11]. Second, S p at the 
heavily-doped n + surface is likely dominated by temperature insensitive Auger 
effects [11]. Third, results similar to those shown in Fig. 2 were also 
obtained from the devices for which passivation by SiO ? sharply decreased S . 

c p 

The fact that I<jq is not negligible below ~ 100 K indicates that most of 
the minority holes are in the extended band states at T ~ 100 K with only a 
small fraction trapped at the tail states and immobilized. For the 
excitations used to measure I^T) ant * also f° r low-level hole injection in 
the dark used in [1], the hole concentration in the n + region is very small. 

We can then conjecture that the penetration of the localized tail states into 
the energy gap is only a very small fraction of the energy gap. This 
supposition, based on our data, agrees with recent theoretical calculations of 
the band structure of heavily doped Si [12]. 

It follows directly from our trapping model that p-(minority) < 
jji(majority) could apply also in the absence of the tail states because shallow 
impurity levels can act as traps near the minority-band edge (Fig. 1). This 
suggests that uOninority) may strongly depend on compensation. A lack of 
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compensation in their epitaxial p/n diodes may be responsible for the 
experimental observation by Dziewior and Silber [13] that p(minority) = 
p(majority) for concentration in Si < 10 19 cm -3 . 

Although the assumptions needed to derive (7) from (5) introduce some 
inaccuracy in the model for the activation energy of minority-carrier mobility 
(diffusivity), the general conclusions derived by comparing (7) with the 
experimental trends exhibited in Fig. 2 are anticipated to remain valid. The 
exact value of the activation energy in a single-trap model remains in 
question. When more detailed knowledge concerning the values and the 
temperature dependencies of parameters in (5) become available, one can then 
use data presented here to explore aspects of the minority-carrier band tail. 

The purpose here is less ambitious. We have demonstrated activation 
behavior of the minority-carrier mobility and diffusivity, exhibiting thereby 
an activation energy of the order of kT (for 30U K) as one expects for the 
prominently active tail states. Moreover, our data in Fig. 2, as interpreted 
here, supports the inadequacy of the commonly used assumption of equal 
majority- and-mi nority-carrier mobilities; in turn, this supports 
interpretations, such as that in [1], of larger values of energy-gap narrowing 
than is common in the literature. 

Lastly, we point out that the assumed value, m = 0, and values m < 0 in 
(4) are self-consistent with the experimental results of Fig. 2, but values 
m > 0 are not. However, regardless of the value of m, the I sc T m vs * plot 
shows an activation behavior, thus the assumption of p(majority) = p(minority) 
is always wrong (see Eq. 8). 
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FIGURE CAPTIONS 


Fig. 1. Qualitative illustration of the band edges of heavily doped 
n + -silicon. The broken lines show the unperturbed parabolic 
bands. The positions of both the electron and hole quasi-Fermi 
levels is also indicated. The arrows near Ey indicate hole capture 
and emission by the tail states and by the acceptor level from the 
p-type substrate. The penetration of the tail states into the 
forbidden gap is assumed to be very small in comparison with the 
bandgap E Q = E^ - Ey. 

Fig. 2. Normalized short-ci rcuit photocurrent versus 1/T for \ = 0.4 pm. 

The illumination density for curve 1 is five times larger compared 
to that of curve 2. 
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CHAPTER SIX 

DETAILED INTERPRETATION OF PHOTOCONDUCTIVITY 
DECAY RESPONSE FOR LIFETIME DETERMINATION 

Our purpose here is to clarify the relation between the observed 
transient photoconductive decay and a previous development of the detailed 
physics [1] underlying this transient. 

To make this clarification concrete, we begin with a review of the work 
of Othmer and co-workers, who compared photoconductivity-decay lifetime with 
steady-state lifetime for silicon [2]. Figure 1 reproduces this comparison. 
The steady-state lifetime was measured by steady-state photoconductivity and 
the transient lifetime was measured by transient photoconductivity. By 
introducing various impurities (Zr, Mn, Mg, Ni , and V) known to produce 
recombination centers in Si and by varying the concentrations of these 
impurities, Othmer and co-workers varied the steady-state lifetime. Note the 
scatter of the data. Note also the large number of data points yielding a 
photoconductivity decay lifetime considerably more than a factor of two 
greater than the steady-state lifetime. 

From Othmer's results we see that the photoconductivity decay is 
characterized by an apparent single characteristic time; yet this 
characteristic time differs considerably from the steady-state lifetime. 

To suggest an explanation for this apparent anomaly, we review briefly 
some aspects of our prevous theoretical work on fast photoconductive switches, 
which we have called photoconductive circuit elements [1]. In these devices, 
incident light creates hole-electron pairs in very high injection in a nearly 
intrinsic region of a p/i/n or n/i/n or p/i/p photoconductor. Though a 
combination of a Dember field and a high surface recombination velocity causes 
a peak in the hole and electron density in the y direction normal to the 
surface on which the light shines [3], integration of the continuity equation 
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in the y direction leads one to analyze a related device having no dependence 
of variables in the y direction in the active region. (This approach 
resembles that used in MOSFET theory to reduce the number of space variables 
required.) 

Then, provided the injection level is high enough that the particle 
dielectric relaxation time is much less than the drift transit time in the 
x direction resulting from an applied electric field, the device is in a 
quasi-neutral condition [1]. Thus the turn-off transient becomes a problem in 
quasi-neutral photoconducti ve decay in the simplest model: that in which the 
x or y coordinates do not enter. 

Despite this reduction to a simple model, some complexities arise from 
recombination through many traps--a prominent feature of this class of 
photoconducti ve device [1]. For example, for incremental variations in the 
particle electrochemical potentials, variations that are much smaller than a 
thermal voltage, a detailed description of decay involves two characteristic 
times for each participating bound state in the forbidden band. This 
contrasts with the single characteristic time, the high-level recombination 
lifetime, that prevails for the time-invariant steady state. Note that these 
multiple characteristic times apply only for incremental times. In Ref. 1, 
apparently the first method is developed for determining large-signal 
transient that incorporates these multiple time constants. 

For the recombination parameters indicated in Fig. 2, we compare the 1/e 
fall time calculated by this method with that resulting from the standard 
steady-state lifetime, a single-level Shockley-Read-Hall mechanism. From this 
comparison we see two key results. First, the single-lifetime calculation 
underestimates the 1/e fall time by about a factor of 1.5 in this example. 

This effect results from the trap readjustment time [1], which is neglected in 
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the standard single-lifetime model. Second, the net effect of two 
characteristic times, which exist for incremental times, yields an essentially 
single-exponential decay for the large-signal response. 

The combination of these two results suggest an interpretation of the 
data of Fig. 1 different from that advanced by Othmer and co-workers [2]. We 
suggest that the observation of photoconduct ive decay exhibiting an apparent 
single characteristic time that differs considerably from the steady-state 
lifetime is not an anomaly. Rather this results directly from the underlying 
physics of photoconductive decay for samples having a large density of 
recombination centers. 

This physics has several features. First, for incremental times and 
small -signal responses, as defined above, two characteristic times are needed 
for each recombination-center level to describe fully the transient decay. 

This detail of the physics was first advanced by Sandiford [4] and later was 
extended by Sah [5]. Second, the incremental decay response can be used to 
yield the large-signal decay response. A method for doing this appears in 
Ref. 1. Third, the large-signal response has an apparent single 
characteristic time (which exceeds the steady-state lifetime) that prevails 
over several decades of decay of the excess carrier densities. This third 
feature was overlooked in Ref. 1, and apparently has not appeared previously. 
It provides a key to the detailed interpretation of photoconductive-decay 
experiments. 

We acknowledge the support of the Silicon Solar Array Project of the Jet 
Propulsion Laboratory, California Institute of Technology, in conjunction with 
NASA and DOE. We thank Bob Hammond and C. Tang Sah for help on technical 
aspects of this and related work. 
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FIGURE CAPTIONS 


Comparison of photoconductivity decay lifetime measured by 
Westinghouse and steady-state lifetime measured by Northrop on 
symbols derived from the same ingots [2]. 

Comparison of decay times for capture rates (thermal velocity times 
capture cross section) c n = c p = 10" 8 cm 3 /s, Nyy = 10 18 cm 3 . 

Curve (a), calculated by using two transient time constants, 
contrasts with curve (b), calculated using a model based on a 
single steady-state lifetime. The 1/e decay time is 1.5 times 
longer than that predicted by a single steady-state lifetime for 
the parameters chosen. This theoretical dependence, for InP, will 
scale properly, for example, for Si, where a capture cross section 
of 10-15 cm -3 and a thermal velocity of 10 8 cm/s (> 10^ cm/s for 
electrons in InP) and a flaw density of lol 3 will yield % - 1 
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CHAPTER SEVEN 

GENERALIZED RECIPROCITY THEOREM FOR SEMICONDUCTOR DEVICES 

In a recent paper, Donolato 1 has proven a version of the reciprocity 
theorem for the charge collection by a p-n junction in the presence of a unit 
point generation source. The proof was restricted to a uniform material in 
which only the lifetime was allowed to vary with position. We extend this 
proof to p-n junction devices having three-dimensional space dependencies of 
all material parameters. In addition, we prove that reciprocity holds also 
for transients and state another form of reciprocity relating quasi-Fermi 
levels and point generation in a quasi neutral region under low injection 
conditions. 

The conventional reciprocity theorem relating input voltages and output 

2 

currents in the base of a transistor was proven by W. Shockley et al . in one 
of the early papers on p-n junction transistors. In this letter we generalize 
Shockley's procedure in order to prove reciprocity for charge collection. 

The principal quasineutral region of the diode, where carrier generation 
and recombination takes place, is bounded by the junction Sj and the rest of 
the surface S L (lateral and back surface) characterized by a position 
dependent surface recombination velocity s'. The region is assumed to be 
p-type in low injection. Extension to n-type regions is straightforward. Two 
different excitation cases are considered. In case 1 the junction is forward 
biased by a voltage V in the dark. In case 2 the junction is short circuited 
and a point generation source of strength G is applied at point r' . If J, n 
and J' ,n' are the minority (electron) current density and excess minority 
carrier concentration for cases 1 and 2 respectively, then the boundary 
conditions for each case are 

n((“-]) q\J 

rpry = exp( *r> * i ■ n '‘ r j> ■ 0 • n> 

0 J 
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where pj is a point on the junction edge Sj and n 0 is the equilibrium minority 
carrier concentration. The continuity equation for cases 1 and 2 is 


1 

q 


v*o 



- — + G 6(r - i*' ) 

X v 


( 2 ) 


_ where 6(r - r‘ ) represents a unit point carrier generation source applied at 
r = r' , and i is the lifetime. Because of low injection, the electric field 
E, the mobility \x and the lifetime t are independent of the excitation. The 
system is linear. 

Let vector A be defined by the equation 


A = S- O' - 5- 0 . 

"o "o 


The surface integral J A»ds vanishes since on S|_ 


A*dS = -q(2$— s' .^s'l =0 . 


n n 


From Eq. (1) this integral on Sj becomes 


(3) 


' (4) 


J A*ds = [exp(-^y) - lj J J 1 *ds 


[exp(-jjy) 


i] r 


(5) 


where I* is the current due to the point source exiting the p-type region 
through the junction. Therefore, if S is the entire surface enclosing the 
quasineutral region, then from Eqs. (A), (5) the total surface integral of A 
becomes 


J A»ds = J A*ds + J A*ds = [exp(^L - 1] I' . (6) 

S S L S J 
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By using Eq. (2), the Einstein relationship for non-degenerate minority 
carriers, D = (kT/q)p, and the expression for the electric field E = 
-(kT/q)vln n Q , the divergence of A can be written as 


V»A = v*C~ J'] - V-C~ J3 ~ ~ Cvn - nvln n ] • J ' + 
n o 1 o n o 

jj- V * J ' " rT" £ vn ‘ " n ' v1n n o> J ' n~ V ’ J = 
oo n o 


—■ [Vn + nE]*[qpn' E + qDvn' j 


q ~ + G 5(r 
n o T 


r' ) J 


-jf- [Vn* + n'EMqpnE + qDvnj + q [- -J] = 
o o 


i-nkTCvn tJ T .E].[l f » l E + vn'] +q ^ - q j-G 6(r - r') 

0 0 0 

- ^ pkTCvn 1 +^r n'E].^ nE + Vn] - q 0 = - q G S(r - r‘ ) . (7) 


By using Guass' theorem and from Eqs. (6) and (7) 


J A*ds = / V»A dv => [exp(^-) - 1]I' = 
S V Kl 


- q / 7 T S 6(r - r 1 )dv 
V o 


n(r') 

n JFJ 

o' 


(B) 
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Equation (8) is similar to Uonolato's final expression. It is more general 

because we include position dependencies of p., 0 and n Q (x). 

We introduce the minority carrier quasi-Fermi level in dark, F n , 

referenced to the majority carrier Fermi level, which is assumed position 
3 

independent . Then, Eq. (8) becomes 

exp(F (r‘ )/kT) - 1 

V - - exp(qV/kT) - 1 G * ^ 

because of linearity. Thus we can choose V conveniently so that 
exp(qV/kT) -1=1. Then Eq. (9) becomes 


I' = - q [exp(F n (r')/kT) - 1] G . (10) 

Equation (10) expresses that the probability of collection of a minority 
carrier generated at point r* is exp(F n (r' )/kT) - 1. 

In a similar manner, one could prove another form of the reciprocity 
theorem applying to a quasineutral region in low injection: the minority- 
carrier quasi -Fermi level at point 1 due to a point source applied at point 2 
is the same as the minority-carrier quasi -Fermi level at point 2 due to the 
same source applied at point 1. 

These theorems hold also for transients. This is proved by replacing n,I 
and n',r by their Laplace transforms and by transforming the continuity 
equations for zero initial conditions: 


dn 

at 



^ => sn(s) 


\ V.J(s) 


njs) => 
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(11a) 


- i V-J(s) = - (s + ^)n(s) , 


= I vJ' - — + G 6(r - r' ) => 
5t q t ' 


- i v-o' (s) - - (S + i)n(s) ♦ S(s)6(r - r' ) (lib) 

Substitution of Eq. (11) into (7) and (8) will give the same result as in the 
steady state. 

Equation (10) could be useful in solar cell theory since it allows the 
calculation of short-circuit current from the analysis of the device under 
dark conditions: 

Ic r = l (- q) / Cexp(F (r)/kT) - 1] G(r)d 3 r , (12) 

bt i=1>2 i m 

where 1,2 represent the quasi neutral emitter and base of the cell, F m i s the 
minority carrier quasi-Fermi level in each region in dark and under bias V 
such that exp(qV/kT) - 1 - 1 and G is the generation rate. Equation (12) 
holds under the same assumptions underlying the proof of the reciprocity 
theorem and ignores generation in the space charge region. 

The major assumption in deriving Eq. (8) was low level injection so that 
the system could be treated as linear. Low injection guarantees linearity of 
the system even in the presence of the effects of heavy doping or graded 
chemical composition. In this case, in the equation -vln n Q = qE/kT, E will 
be the minority carrier field which includes the quasi-fields arising from the 
position dependent energy gap electron affinity and density of states. 
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Therefore, the original Shockley's proof of reciprocity can be extended 
to heavily doped three dimensional base regions to obtain relations between 
input voltages and output currents. In this sense, Shockley's proof is more 
general than the recently published proof of reciprocity in one- dimensional 
transistor base regions which takes into account heavy doping effects. 

The practical significance of reciprocity is that it can reduce the 
number of unknowns in a problem and reduce the amount of calculations 

3 

needed. As an example, in the recent paper by del Alamo and Swanson , the 
quantum efficiency of an emitter could have been found by integrating 
Eq. (10), since the minority carrier Fermi level as a function of position was 
already numerically calculated under dark conditions. In addition, their 
forward and reverse transport factors and saturation currents are related by 
the relation = a r J op . Thus the present reciprocity theorem can find use 

in the measurement of such device parameters as the surface recombination 
velocity at a polycrystalline silicon contact. 

Although not discussed here, reciprocity relations enter the physics of 
many different phenomena, such as heat conduction or thermoelectric effect. 
These relations are based on fundamental considerations relating to detailed 
balance under thermal equilibrium and Onsager's principle of microscopic 
reversibility. 
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CHAPTER EIGHT 

MINORITY-CARRIER DIFFUSION COEFFICIENTS 
AND MOBILITIES IN SILICON 

Minority-carrier mobility is one of the principal parameters determining 
the transport properties of minority carriers in p-n junction devices. The 
majority-carrier mobility (conductivity mobility) is important for 
determination of the resistance of semiconductor regions. Accurate data and 
theoretical expressions for the conductivity mobility in both the n- and p- 
type Si are available [ 1 - 6 ]. Only very limited date are presently available 
for the minority carriers. Prince [ 7 ] measured the minority-carrier drift 
mobilities of electrons p n and holes Up using the Haynes-Shockley pulse 
method. Dannhauser [ 8 ] and Krausse [ 9 ] measured the sum (p n + pp) from an 
ohmic drop of forward-biased p-i-n diodes. Recently, Dziewior and Silber [ 10 ] 
determined the minority-carrier diffusion coefficients, D = (kT/q)p ,in the 
doping range from about 10^ cm“^ to about 10^ cm“^ by measuring the complex 
diffusion length of the minority carriers generated by 10.7 MHz optical . 
excitation. These recent data show good agreement with the conductivity 
electron mobility, but the hole mobility in n-type Si was up to ~ 30 % higher 
than the hole conductivity mobility in p-type Si. Burk and de la Torre [ 11 ] 
confirmed the result of [10] for one value of doping (Npp = 2.6 x 10^ 
cm"3). In contrast to the results of [ 10 ], the measurements of Prince [ 7 ] 
show that both p n and Pp are smaller than the majority carrier values. Some 
possible reasons for the discrepancies are discussed in [10]. Additional 
reasons may be due to the insufficient purity of Si crystals used in 1954 by 
Prince [ 7 ], 

Because of the scarcity of the accurate and reliable values for minority- 
carrier mobilities and because of discrepancies among the results quoted 
above, modeling of p-n junction devices is based on a convenient assumption of 
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equal minority and majority carrier mobilities. In addition to [10], several 
other theoretical [12-14] and experimental [15,16] works suggested that the 
above assumption is not valid, particularly in heavily doped Si. To verify 
this assumption, we present a new and accurate method for measurement of the 
minority-carrier diffusion coefficients in semiconductors. The method is 
applicable for both low and highly doped materials and is demonstrated here 
for low doped Si . 

The present method is based on a direct measurement of the transit time 
of the minority carriers through a narrow region of a p-n diode. Assume a p + - 
n diode in which the hole diffusion length L p W, where W is the width of the 
n-base region terminated by an ohmic contact with recombination velocity S. 

As was shown by Gonzalez and Neugroschel [17], the real part of the small - 
signal hole admittance for frequencies to ^ 10/x p follows the w 1 / 2 dependence 
according to 

6 p F = K(y L p )(<V 2)1/2 (1) 

o 

where K = (q/kT) ( qA n.j/N nn ) exp [ (qV/kT ) -1] is a constant depending on device 
area A, base doping density Nqq, temperature n^ 2 (T), and applied voltage V, 
and r p = (L p /n p ) is the hole lifetime. The extrapolation of this dependence 
to the low-frequency value for large S p given by the standard relation 
G p F = K(0p/L p ) coth(W/L p ) gives an intercept at 

oj = (2D p /Lp ) coth 2 (W/Lp) . (2) 

This can be further simplified for (W/L p ) ~ 0.3 as 
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(3) 


“i - 2 V« 2 = 1/T t 

where x. = W/2D is the transit time across the base, 
t p 

Equations (2) and (3) are the basis of the method proposed. The method 
requires the measurement of only one parameter gJ^, and is applicable for any 
narrow region with known value W, regardless of doping. 

To demonstrate the method, we show in Fig. 1 the measured and theoretical 
dependencies of Gp vs frequency f = w/2ir for a p + -n device with Nqq = 7.25 x 
10^ cm"3 and W = 310 ym. The intercept frequency is fj = 6.59 kHz. To 
obtain D p we use (2) with accurately measured [17-19] L p = 500 - 10 ym and 
obtain D p = 15.71 - 1 cm^/s and y p = 608 - 40 cm^/Vs at T = 300 K. This value 
for the mobility is about 30% higher than the majority-carrier value of 470 
cm^/Vs [2], The error of about ± 40 cm^/Vs is an estimated total experimental 
accuracy which is discussed in more detail below. 

The error in using (2) instead of (3) is very small. For example, a 10% 

error in L p for the sample in Fig. 1 yields less than 1% error in D p . The 
error can be eliminated by using thinner samples as shown later. In obtaining 
the above result from the measured conductance we assumed a negligible 
contribution of the p + -emitter to the conductance. This can be justified by 
noting that the intersect frequency corresponding to the emitter transit time 
t tE is f IE = (1/2tt) ( l/t tE ) = (1/2-n) (2D n /w|) = 10 9 Hz for W E = 0.2 ym and D n ~ 

1 cm^/s, which is much smaller than the observed base f IB = 6 kHz in Fig. 1. 

Note also a small hump at the knee in Fig. 1, which is characteristic of the 
base* and would disappear if the emitter (or shunt conductance) contributions 
were larger than about 5%. The dominance of the base was in addition 
confirmed by comparing the measured dc current with that calculated from the 
base parameters. For f ~ 60 kHz, the G vs f dependence for the sample in 


<c - ^ 
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Fig. 1 begins to bend upwards and converges to the G « dependence because 
of the effects of the series resistance of the sample. This range is not 
displayed in Fig. 1. 

In the samples where the emitter contribution is small, but not 
negligible, where m < 2 at the knee, but for higher values of to 

(where the base contribution prevails) it converges to G^ « co^2. We can 

then use (1) in the tu ^ ran g e to obtain 

n p = (1/irf) (G HF /K) 2 . (4) 

Using (4) for the device from Fig. 1 we obtain y p - 677 cm^/Vs which is in 

i o 

very good agreement with the 608 - 40 cm £ /Vs value obtained from (2). The 
latter value is more accurate because of large sensitivity of K in (4) to 
n i (Dp * 1/K * 1/rK ). The value of used above was from Putley and 

Mitchell [20]: n? = (9.61 x lO^jT^exp [-1206/(kT/q)3. Note also that (4) is 
valid regardless of the value of S at the contact. 

Table I summarizes the results from the devices used in this paper. The 
conductance was measured using a Wayne-Kerr B224 bridge and Hewlett-Packard 
4274A and 4275A LCR meters. The device BJT is a bipolar n + -p-n + transistor 
with a uniformly doped base with VI = 16.5 ym and L n = 115 ym [21]. The value 
for y R for this transistor was obtained from (3) by measuring the frequency 
dependence of i c /v BE , which is equivalent to (1) for diode. The advantage of 
using a bipolar transistor is that the collector current i c always follows the 
ideal 'exp (qV BE /kT) dependence and that the i c depends on the base properties 
only. The base width W was measured using a precision mechanical indicator 
with a resolution of ~ 0.5 ym on the diode structures and by angle lapping and 
staining for the transistor. The error in W was less than about 3% for diodes 
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(W = 240-310 um) and was due mostly to the nonuniformity of the large area 
samples used (~ 0.5-1 cm 2 ). This error can be reduced using smaller size 
samples. The dopant density was obtained by the junction C-V method. The 
estimated overall error of the method for the low-doped devices showing a good 
G « f^2 dependence is about 5%. The error will be larger for thinner high 
doped devices due to the increased error in determining W. 

Figures 2 and 3 show the results from Table I plotted together with the 
results from [10] and one point from [11]. For comparison we also show the 
majority-carrier values by Thurber, et al . [1,2]. Our present data show that 
the electron minority-carrier mobility is about the same as the majority- 
carrier mobility, but the minority-carrier hole mobility is larger than the 
corresponding majority-carrier mobility. For N Q0 ~ 10 15 cm” 3 the difference 
is about 30%. Our data are thus consistent with the previous work of Dziewior 
and Silber [10]. Temperature dependence of y p and y n was also studied on 
samples from Table I in the range 300K-360K. The preliminary results are that 
y * T" m , where m - 2 - 2.5, similar to the temperature dependence of the 
majority- carrier mobilities [5,6], 

In summary, a new method for direct measurement of the minority-carrier 
diffusivity and mobility was presented and demonstrated for low-doped Si. The 
method is accurate and requires only a knowledge of geometry in properly 
prepared test samples. Significant differences (~ 30%) between the minority 
and commonly used majority carrier values were obtained for y p , while y n was 
found in good agreement with the majority carrier values. 

The measurements 'can be extended to higher concentrations up to 
~ 10 20 cm” 3 . This however requires measurements at higher frequencies, 
because W has decrease as concentration increases to assure L ~ W. Results 
for the concentration range 10 19 - 10 20 cm” 3 will be reported separately. 
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TABLE I 


Summary of results at 300 K. 


Device 

Base 

Dopant 

Base 

Doping 

(cm- 3 ) 

y P 

(cm 2 /Vs) 

y n 

(cm 2 /Vs) 

S2 

Phosphorus 

7.25 x 10 14 

608 ± 40 


5N 

Phosphorus 

3.32 x 10 15 

620 ± 40 


SP9 

Boron 

1.34 x 10 15 


1288 ± 50 

BJT 

Boron 

1.40 x 10 15 


1256 ± 50 

5P 

Boron 

6.31 x 10 15 


1200 + 100 
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FIGURE CAPTIONS 


Normalized small-signal conductance versus frequency for the sample 
S2 from Table I. The conductance follows the G * f^2 dependence 
for f ~ 25 kHZ. 

Minority-and-majority-carrier hole mobilites versus doping 
concentration at 300 K. The error bars for the present data are also 
shown. The solid line is a computer curve fit of the majority- 
carrier data given by Thurber et al . [2], The broken line is an 
approximate fit to the present and previous minority-carrier results 
[ 10 ]. 

Minority-and-majority-carrier electron mobilities versus doping 
concentration at 300 K including error bars for the present data. A 
computer curve fit for majority electrons by Thurber et al . [1] is 
shown for comparison. 
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CHAPTER NINE 

DISCUSSION AND RECOMMENDATION 

Heavily doped polysilicon (Chapters Two and Three) provides the vehicle 
for a new class of silicon solar cell designs. Ultimate efficiencies may 
exceed 23%. This view is consistent with the projections of C. T. Sah (see 
Chapter Two) made in his research supported by JPL. 

Our study described in Chapter Four strongly suggests the validity of the 
Auger-impact rate coefficients for Auger band-band recombination determined by 
Dziewior and Schmid (1977), despite the views of other workers. Chapter Five 
provides support for the model for low minority-carrier mobilities in highly 
doped silicon that was proposed by Neugroschel and Lindholm in 1983, under the 
support of JPL. 

All of the issues described in Chapters Two through Five are critical 
both to the design of high efficiency silicon solar cells, with the aid of 
computer programs (such as JPL's SEEMA programs), and to the assessment of the 
maximum obtainable practical conversion efficiency for silicon solar cells. 

Chapter Six may help prevent a misinterpretation of the connection 
between recombination lifetime of starting material and the photoconducti ve- 
decay data by which it is measured. 

The generalized reciprocity theorem of Chapter Seven is presented in 
brief mathematical style. Its applications are wide, and we anticipate that 
many useful design results will flow from it. The connection it provides 
between dark current and internal quantum efficiency is outlined as a prelude 
to further uses. 

Minority-carrier mobility and diffusivity in dilutely doped silicon has 
received little attention since the work of Prince in the 1950's. Appropriate 
values have been assumed to be available in the standard tables found in text 
and reference books. Our new results should be incorporated into the 
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numerical values used in either analytical or computer-aided design of silicon 
solar cell . 

More important than the results of Chapter Eight is the method by which 
the results are obtained. An extension of this method can yield the minority- 
carrier diffusion coefficient and mobility of heavily doped silicon. This 
determination is essential for systematic design of silicon solar cells. It 
will lead also to a determination of the surface recombination velocity 
adjoining highly doped layers. The work of the Stanford group (del Alamo, 
Swanson, and Swirhun) supported by OPL has clarified the utility of * 

considering the minority-carrier equilibrium concentration and the effective 
doping concentration as an inseparable, and experimentally determinable, 
quotient. This utility extends only to devices for which the heavily doped 
layers end at surfaces having very large surface recombination velocity S. 

When the surface is passivated, the minority-carrier diffusivity must be 
determined to yield S, and S and effective doping concentration must 
apparently be separately determined experimentally. Because most solar cells 
will have such passivation, this becomes a key issue. Furthering of the 
method of Chapter Eight will be part of our approach toward the solution of 
this design and analysis problem. 
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FIG. 3 
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